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ABSTRACT 


— — ^  A  chemically  reacting  shear  layer  between  various 

gases  was  investigated  in  a  new  type  of  blow-down  wind 

tunnel.  The  gas  streaas  were  inert  (helium,  nitrogen  or 

ergon),  hut  carried  up  to  10Z  concentration  of  reactants, 

one  being  osone  and  ^e  other  nitric  oxide.  The  resulting 

+  NO  NO^7  +  <T,  was  essentially  diffusion 
2  2 

limited  and  spontaneous,  enabling  the  temperature  rise  to 

be  varied  at  will  from  sero  up  to  20d®C  mean,  Flows  of 

I db  I?'*  t 

Reynolds  number  up  to  5  x  (o**  were  investigated.^^ 

Simultaneous  shadowgraphs  at  two  different  wavelengths 
(blue  and  middle  ultraviolet)  were  taken,  the  blue  light 
shoving  the  product  NO^  and  the  UV  showing  the  reactant  0^. 

Mean  velocity  and  temperature  profiles  of  the  layer  were 
obtained  for  various  reactant  concentration  ratios  and 
temperature  rises. 

The  effect  of  heat  release  by  the  chemical  reaction 
was  to  strongly  reduce  the  entrainment  rate  of  free  stream 
fluid  into  the  shear  layer,  and  to  apparently  inhibit  the 
smallest,  high  frequency  turbulence  scales.  The  effect  of 
changing  the  reactant  concentration  ratio  was  to  skew  the 
mean  temperature  profile  towards  the  minority  reactant  side 
of  the  layer,  but  the  dependence  was  not  strong. 

The  measured  temperature  profiles  were  also  compared 
with  those  ealculsted  from  probability  density  functions  (PDFs) 
for  the  mixture  composition  at  various  points  across  the 
layer,  which  were  measured  in  a  non- reacting  turbulent  shear 


reaction. 


P 

3 


layer  by  Konrad  (1977).  Good  agreement  was  found  within 
the  temperature  rise  range  of  the  experiment,  especially 
for  tenperature  rises  lass  than  100°C.  A  procedure  for 
inverting  the  above  calculation  (i.e. ,  finding  the  PDFs 
from  measurements  of  mean  temperature  rise  profiles)  was 
used  to  produce  a  set  of  PDFs  comparable  to  Konrad's  which 
would  predict  temperature  rise  more  accurately. 

The  strong  effect  on  the  temperature  rise  of  changing 
the  density  ratio  of  the  two  free  streams  was  demonstrated 
and  shown  to  be  in  good  agreement  with  the  effect  predicted 
from  Konrad's  non-reacting  PDF  measurements.  It  was  also 
contrasted  with  the  very  different  result  obtained  from  an 
eddy  dlffuslvlty  model  of  the  turbulent  shear  layer. 

Experiments  with  an  argon  shear  layer  which  consistently 
produced  unexpectedly  low  temperature  rises  are  reported. 
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1.  INTRODUCTION 


1. 1  Background 

Modern  theory  of  the  fluid  mechanics  of  combustion 
really  began  in  1928  with  Burke  and  Schumann.  Prior  to  that 
time,  there  appears  to  be  no  serious  attempt  to  bridge  the 
gap  between  the  chemists,  with  their  Lagrangian  (time  based) 
idealizations  of  molecular  reaction  rates,  and  the  kinetlclsts 
with  their  practical  reacting  flow  problem,  often  in  the 
Eulerlan  (special)  framework.  Working  with  non-premixed 
combustion  in  laminar  flow,  Burke  and  Schuaann  realised  that 
the  chemical  reaction  between  two  fluids  was  confined  to  a 
zone  at  the  Interface  which  they  called  a  "flame  sheet". 

Reactants  were  fed  to  each  side  of  the  flame  sheet  by  Fick's 
Law  diffusion.  Thus  they  Introduced  the  idea  of  a  "diffusion  - 
limited"  reaction  rate.  Damkohler  (1936)  then  drew  attention 
to  the  influence  of  the  ratio  a,  between  a  characteristic 
diffusion  time  and  a  characteristic  molecular  chemical  reaction 
time.  Where  this  is  zero,  the  flow  is  non-reacting  (frozen 
chemistry)  and  where  It  approaches  Infinity,  the  reaction 
zone  becomes  vanishingly  thin  as  the  reactant  molecules 
cannot  co-exist  in  the  same  region.  This  second  case  represented 
the  diffusion  limited  flame  sheet  of  Burke  and  Schumann,  and 
has  received  the  greatest  attention  in  the  literature  ever 
since  because  the  chemical  reaction  rate  term  can  be  eliminated 
from  the  flow  equations,  which  results  in  a  considerable 
simplification  over  any  intermediate  value  of  a.  Friedlander 


and  Keller  (1936)  analytically  solved  the  static  isothermal 
plana  diffusion  layer  with  a  one  step  second  order  reaction. 
The  theory  of  combustion  in  laminar  flows  was  reviewed  by 
Williams  (1971). 

Application  of  laminar  flow  concepts  to  turbulent  flow 
increased  markedly  in  the  1940's.  Hawthorne.  Hottel  and 
Weddel  (1949),  working  in  a  jet,  measured  the  mixture 
fraction,  c,  the  molar  fraction  of  fluid  originating  from 
the  jet  in  any  sample.  They  showed  that  the  product  con¬ 
centration  could  not  be  related  to  the  time  average  mixture 
fraction  at  a  point.  What  was  required  was  a  complete 
probability  density  function  for  the  mixture  fraction  at 
that  point.  (This  idea  is  analytically  described  in 
Appendix  A).  Their  work  underlined  the  Importance  of  the 
concept  of  molecular  scale  mixing  and  unmixedness  in 
diffusion  limited  turbulent  flows.  The  time  average  turb¬ 
ulent  system  differed  from  the  non-turbulent  in  the  fact 
that  the  Interface  could  not  be  a  thin  surface,  even  with 
Infinite  reaction  speed.  The  necessity  for  a  detailed 
knowledge  and  understanding  of  the  turbulent  fluid  mechanics 
was  thus  confirmed,  but  much  significant  work  at  that  time 
occurred  in  the  field  of  homogeneous  turbulence.  Since 
Kolmogorov  (1941),  it  was  well  known  that  the  small  scales 
of  turbulence  were  isotropic  and  independent  of  the  large 
scales  and  turbulent  energy  was  cascaded  from  the  large 
eddies  to  the  small,  dissipating  eddies.  Hence  the  study 
of  the  homogeneous  turbulence  situation  had  a  role  to  play 
in  the  field  of  shear  flow  after  breakdown.  Some  of  the 


workers  in  this  ares  wars  Batchelor  (1952)  on  the  stretching 
of  fluid  lines  end  surf aces  In  isotropic  turbulence,  O'Brien 
(1971),  and  Corrsln  (1958)  on  the  statistical  behaviour  of 
a  reacting  mixture  under  such  flow  conditions. 

Prediction  of  reaction  ratee  In  turbulent  flow  wee 
approached  by  Toor  (1962),  who  used  an  analogy  between  a 
scalar  species  in  a  non-reacting  flow  and  a  reaction  invarlate 
species  In  the  reacting  flow.  With  the  assumption  of  dilute 
reactants,  equal  dlffusivltles,  and  a  rapid.  Irreversible 
reaction,  the  tine  average  behaviour  of  the  reacting  system 
could  be  deduced  If  the  RMS  concentration  fluctuations  and 
time  average  behaviour  were  known  in  the  non-reacting  system. 
Lin  and  O'Brien  (1974)  extended  Toor's  analysis  to  shear  flows, 
solving  for  the  reacting  species  concentrations  In  a  jet  using 
data  from  Tutu  and  Chevray  (1973)  of  temperature  measurements 
In  a  heated  non-reacting  jet. 

In  the  meantime,  the  behaviour  of  non-reacting  turbulent 
shear  flow  structure  had  received  enormous  attention.  The 
first  detailed  study  of  the  shear  layer,  using  hot  wire 
anemometry.  Is  attributed  to  Liepmann  and  Laufer  (1947). 
Vfygnanskl  and  Fiedler  (1970)  repeated  the  study.  Including 
higher  order  velocity  fluctuation  products  and  intermittancy. 
Some  pressure  data  was  added  by  Spencer  and  Jones  (1971) 
using  a  special  pressure  transducer.  The  recent  recognition 
of  large  structure  in  the  shear  layer  (Brown  and  Roshko,  1974; 
Winant  and  Browand,  1974;  Dlmotakls  and  Brown,  1976)  has  led 
to  a  much  deeper  understanding  of  the  entrainment  and  mixing 
processes.  The  important  experiments  of  Konrad  (1977)  who 
used  an  aspirating  probe  (Brown  and  Rebollo,  1972)  to  measure 
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species  end  density  elucidated  the  highly  non  gaussian  PDFs 
of  the  mixture  fraction  In  the  shear  layer  (In  the  absence 
of  better  Information,  model  builders  were,  and  are,  often 
uneasily  forced  to  assume  Gaussian  distributions) . 

At  this  stage  a  number  of  computational  models  exist 
and  appear  to  predict  mixing  and  reacting  processes  In 
turbulent  flows  (Spalding,  1974).  However,  the  models  badly 
need  quantitative  experimental  verification,  and  In  particular 
the  greatest  unknown  la  the  Interaction  between  turbulence 
and  chemistry.  Reacting  experiments  to  date  have  fallen  into 
two  distinct  categories. 

1.  Those  where  conditions  approaching  Toor's  analogy 
apply  vis.  Insignificant  heat  release,  usually 
conducted  specifically  to  Isolate  the  secondary 
effects  of  reaction  (Batt,  1977;  Shea,  1976). 

2.  Those  with  high  or  extreme  heat  release  (causing 
Incandescence)  where  most  of  the  highly  developed 
and  sensitive  instrumentation  of  the  Isothermal 
research  Is  unusable  (Bllger  1976;  Wahl  et  al 
1949  etc.). 

The  present  research  is  the  first  experiment  to  bridge 
these  categories,  having  infinitely  controllable  heat  release 
from  zero  to  a  maximum  of  approximately  200°C  time  averaged. 
1.2  Present  Research 

The  requirement  in  designing  this  experiment  was  to 
abstract  a  flow  which  had  a  vide  applicability  to  real  engin¬ 
eering  situations  (in  particular  the  chemical  laser)  and  at 
the  same  time  be  as  simple  as  possible  for  experimental  and 
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analytic  tractabllicy.  Thera  vaa  little  difficulty  in 
chooeing  the  plane  ahear  layer  (Figure  1.1)  for  the 
study.  Although  it  required  a  dual  flow  facility,  it 
had  the  following  characteristics. 

1.  Two  dimensionality. 

2.  Relatively  fast  achievement  of  a  linear 
similarity  in  the  x  direction. 

3.  Considerable  previous  investigation  (non- 
reacting)  e.g.  Liepmann  and  Laufer  1947, 

Wygnanski  and  Fiedler  1970,  Brown  and  Roshko 
1974,  Konrad  1976. 

A  near  ideal  chemical  reaction  to  use  in  this  study 
was  the  nitric  oxide  -  ozone  reaction. 

NO  +  0  0  +  NO 

3  2  2 

which  has  been  well  documented  (Johnston  and  Crosby,  1951, 
1954;  Herron  and  Hine,  1973),  partly  from  research  motivated 
by  atmospheric  pollution  problems. 

It  is  a  simple,  second  order,  blmolecular  reaction 
accomplished  in  a  single  collision.  Although  other  reactions 
can  occur  in  the  flow,  this  reaction  dominates  within  the 
temperature  range  of  the  experiment.  (See  Appendix  D  for 
details  of  the  chemical  kinetics).  In  particular,  these 
other  reactions  can  occur  in  the  rig 

1.  dissociation  of  NO 

2 

NO*  2N0 
2  4  2 


Inert  Gas  +  $  10%  NO 


Velocity 


Fig.  1.1  Chemically  Reacting  Shear  Laver  Experiment. 
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oxidation  of  >0 


3. 


w  +  »,—  *>, 

dislocation  of  osona 


2/30^0 
3  2 

Data  froa  Varhoak  and  Danlala  (1931)  shows  that  the 

dislocation  of  II  0  to  MO  is  highly  temperature  depend- 
2  4  2 

ent,  and  that  the  reaction  product  at  the  teaperature 

occurring  la  this  Investigation  Is  substantially  NO  ,  not 

2 

M  0  .  With  dilute  reactants,  when  teaperature  rises 
2  4 

are  small,  the  overall  effect  oa  the  flow  due  to  molecule 
number  change  would  be  Insignificant.  Comparison  of  the 
reaction  rate  constants  for  MO  +  0  and  MO  ♦  0  shows 


2  3 

that  the  first  reaction  la  exteaely  alow  and  can  be  Ig¬ 
nored  in  this  flow.*  The  decomposition  of  osone  into 
oxygen  was  minimized  by  appropriate  design,  maintenance 
and  operation  of  the  rig  (see  2.4.)  and  Its  effect  nulled 
by  run  time  aeasursaent  of  the  ozone  concentration  (see 
2.5.4).  Besides  simplicity,  the  ozone,-nltrlc  oxide  reaction 
features  unidirectionality,  spontaneity,  very  high  speed, 
high  exothermoslty  and  conservation  of  molecular  nunber. 
Ignition  and  extinction  problems  are  thereby  avoided,  and 
It  Is  possible  to  Independently  vary  the  heat  release  in 
the  turbulence  froa  zero  to  a  maximum  dictated  by  experimental 
or  physical  constraints. 

Calculations  in  Appendix  D  show  that  the  chemical  reaction 
time  is  fast  enough  for  the  overall  reaction  rate  In  the  shear 


t  Some  oxygen-nitric  oxide  runs  were  made.  Mo  reaction 
was  observed  In  the  test  section,  although  sons  No^ 
was  present  when  the  gas  finally  entered  the  exhaust 

systsm. 
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layer  to  be  considered  diffusion  limited.  A  further  useful 
feature  is  the  possibility  to  independently  measure  or  photo¬ 
graph  either  NO  or  0  due  to  their  separable  and  unique 
2  3 

absorption  spectrums  from  all  other  species  present  (see 
Figures  2.3). 

To  achieve  a  flow  with  the  chemical  reaction  described 
here  required  a  new  type  of  test  facility.  For  safety  and 
economy  gas  quantities  should  be  small,  hence  an  intermittent 
flow  down  rig  was  implicated.  Since  running  times  were  to 
be  kept  small  (seconds) ,  fast  multichannel  data  acquisition 
facilities  were  also  required.  The  unique  system  built  is 
described  in  Section  2.  It  is  highly  versatile  due  to  its 
ability  to  deliver  a  constant  flow  of  any  reasonably  stable 
gas  mixture  compatible  with  Teflon  and  stainless  steel  and 
has  been  used  for  studies  into  non  uniform  density  jets  and 
jet  noise  besides  the  shear  layers  in  this  investigation. 

The  first  results  produced  by  the  rig  were  published  by 
the  author  in  1977. 

1.2.1  Alms  of  the  Research  Prograane 

The  discussion  in  1.1.  Indicates  that  there  is  a  large 
number  of  questions  to  be  answered  and  probably  the  most 
pressing  of  these  relates  to  the  effect  of  the  heat  released 
from  the  chemical  reaction  on  the  turbulent  flow  field.  The 
present  facility  is  Ideally  suited  to  Incrementally  increasing 
the  heat  release  from  zero  by  simply  changing  the  fuel  and 
oxidlser  concentrations,  and  at  the  same  time  effects  of  changing 
the  fuel-oxldlser  ratio  can  be  studied.  These  two  parameters 
appear  later  as  a  temperature  rise  number,  8  *  AT/Ts,  where  at 


la  a  taaparature  rlae  acale  and  T  tha  abaolute  free  stream 
temperature,  and  a  raactant  concentration,  £  "  C0o^c0jjo' 

Tha  most  important  of  tha  affecta  produced  by  the  chemical 
reaction  to  be  atudied  at  thia  stage  of  understanding  of  the 
combustion  problem  are  those  relating  to  the  overall  flow 
geometry,  the  turbulent  structure,  and  the  molecular  scale 
mixing  rate  upon  which  the  reaction  depends.  Hence  the 
experiment  has  three  main  approaches. 

1.  Effects  of  6  and  £  on  the  time  average  growth  rate 
and  entrainment. 

2.  Effects  of  0  and  £  on  the  rate  of  molecular  scale 
mixing. 

3.  Effects  of  8  and  the  internal  flow  structure  of  the 
turbulent  shear  layer. 

Another  important  aspect  of  the  work  is  to  partially 
check  the  contribution  of  Konrad  (1977)  to  combustion  research. 
Konrad's  measurement  of  mixture  fraction  PDFs  in  the  turbulent 
shear  layer  can  be  used,  with  some  assumptions,  to  estimate 
temperature  rises  in  that  flow.  This  experiment  will  provide 

for  the  first  time  a  direct  check  on  the  PDFs  measured  in  the 
non-reacting  flow,  as  well  as  their  applicability  to  the  re¬ 
acting  flow. 
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EXPERIMENTAL  EQUIPMENT 


2. 1  Plow  Production  and  Exhaust  Treatment 

The  Intermittent  blow-down  facility  used  for  the  present 

study  Is  shown  in  Figures  2.1  and  2.2.  The  gas  charge  for 

each  side  of  the  mixing  layer  is  mixed  and  temporarily  stored 

in  a  gas  tight  bag  made  from  Fluorlnated  Ethylene  Polythene 

("FEP  Teflon")  0.125  mm  thick  and  heat  sealed  at  the  joints. 

These  bags  are  contained  within  two  174  litre  steel  vessels 

which  are  pressurized  by  a  much  larger  reservoir  (4.65m3) 

of  compressed  air.  Flow  is  started  and  stopped  by  a  fast 

operating  (300mS)  ball  valve  before  the  reinforced  PTFE  line 

between  the  FEP  bag  and  the  test  section,  and  is  controlled 

by  a  finely  adjustable  multijet  low  noise  metering  valve. 

It  consists  essentially  of  a  movable  piston  in  a  cylinder 

radially  drilled  with  about  250  small  holes  equlspaced  along 

a  helix.  Not  only  does  this  permit  a  closely  linear  variation 

of  area  with  turns  of  the  hand  wheel,  but  it  minimises  the 
>*• 

noise  produced  by  keeping  the  frequencies  high  and  having  a 
large  number  of  discrete  quadrapole  sources.  The  pressure 
ratio  across  this  valve  is  maintained  high  enough  during  the 
run  to  choke  it,  and  since  the  system  pressure  is  maintained 
by  the  large  reservoir  to  within  typically  3%  as  the  bags 
collapse,  the  mass  flux  and  flow  velocity  in  the  test  section 
are  practically  constant.  After  a  short  time  in  the  settling 
chamber,  the  gas  passes  through  the  usual  flow  straightening 
screens  and  contractions  to  the  test  section. 


TO  OTHER  VESSEL 


Fig.  2,1  Blowdown  Facility  (Only  One  Vessel  Shown) 


Fig.  2.2  Blowdown  Rig-  Shear  Layer  Test  Section. 
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The  diffusing  section  of  each  duct  consisted  of  two 
screens  of  SOX  open  area  perforated  sheet,  the  first  bent 
as  shown  in  Figure  2.2  to  approximate  an  equipotential 
plane  where  the  streamlines  are  diverging.  Four  more 
screens,  two  of  30  mesh  and  two  of  100  mesh  effectively 
reduced  the  turbulence.  The  entire  rig  in  Figure  2.2. 
was  made  of  stainless  steel  except  for  the  contraction 
blocks  (polyester  coated  wood),  the  side  plates  (10  ran  per¬ 
spex)  and  the  moveable  side  walls  (glass).  Table  2.1 
summarises  some  specifications  of  a  typical  flow. 

TABLE  2.1 


High  Speed  Side  Low  Speed  Side 


Contraction  ratio 

4:1 

2:1 

Contraction  exit  size 

100  x  25 

100  x  50  mm 

Splitter  plate  included  angle 

3.6° 

1 

Thickness  of  end  of  splitter  plate 

0.02  mm 

Velocity  (nitrogen) 

25 

5  m/s 

Turbulence  level 

0.25 

0.5  XRMS 

Boundary  layer  momentum  thickness  8 
(Thwaite's  method)  at  end  of  splitter 
plate 

0.06 

0.21  mm 

0  u. 

Be  L«  R®  ^ 

100 

70 

Maximum  run  time 

17  sec 

Typical  run  time 

3  sec 

14. 


Included  for  completeness  here  also  Is  a  plot  of  the  velocity 
profile  very  close  to  the  splitter  plate.  (Figure  2.3*).  The 
reason  for  this  Is  that  the  Initial  conditions  of  the  shear 
layer  have  been  shown  to  be  very  significant  in  determining 
Its  subsequent  behaviours.  The  measurements  were  made  by 
traversing  a  constant  temperature  hot  wire  anemometer.  We 
note  that  the  calculated  momentum  thickness  of  the  boundary 
layer  on  the  low  speed  side  Is  In  excellent  agreement  with 
the  velocity  profile  measured  1  mm  from  the  end  of  the  splitter 
plate. 

Any  mixture  of  gases  compatible  with  stainless  steel  and 
Teflon  could  be  used  In  the  rig  and  the  arrangement  supports 
the  safe  and  convenient  use  of  an  unstable  gas  like  ozone 
(see  Section  2.4). 

The  exhaust  gas,  usually  containing  0^,  N^,  NO^  and  0^ 

or  NO,  Is  collected  In  a  large  FEP  bag  at  the  atmospheric 

pressure  for  temporary  storage  while  a  small  pump  scrubbed  it 

through  a  solution  of  sodium  hydroxide.  This  caustic  solution 

catalyses  any  remaining  0  to  0  and  absorbs  and  neutralises 

3  2 

most  of  the  nitrous  oxides. 

The  success  of  the  system  depends  on  the  FEP  bags  staying 
acceptably  gas  tight  and  although  the  bags  could  withstand  re- 
peated  complete  collapses  at  full  system  pressure  (700  kPa) , 
care  had  to  be  taken  not  to  inflate  them  beyond  their  capacity. 
A  special  purpose  FEP  sealing  tool  was  made  to  manufacture 
and  repair  the  bags.  It  consisted  of  two  100  mm  lengths  of 
12.5  mm  diameter  thick  walled  copper  tubing  heated  internally 


VELOCITY  m/s 
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by  electric  elements  and  thermostatically  maintained  at 
about  260°C.  These  were  brazed  to  the  arms  of  a  spring 
loaded  "scissors"  mechanism  and  handgrip.  Tomake  a 
joint,  the  two  thicknesses  of  FEP  were  sandwiched  between 
two  thin  sheets  of  PTFE,  and  then  clamped  for  a  fixed 
time  by  the  sealing  device.  A  little  adjustment  of 
time,  clamping  pressure  and  temperature  soon  produced 
repeatable  joints  which  were  as  strong  in  transverse 
tension  as  the  un jointed  FEP  sheet. 

2 . 2  Data  Acquisition  and  Traverse  Mechanism 

Data  was  acquired  on  a  multichannel  high  speed 
analogue  to  digital  conversion  system  designed  and  built 
within  the  Department  of  Mechanical  Engineering, 

University  of  Adelaide.  The  system  accepts  data  on  up 

i* 

to  16  channels  simultaneously,  4  digital  and  12  analogue.  f 

Each  input  channel  can  be  Independently  sampled  at  a  pre¬ 
programmed  rate  up  to  0.5  MHZ.  The  analogue  inputs  are 
converted  to  digital  values  and  the  numbers  stored  in  a 
shift  register  memory  of  15k  x  8  bit  bytes.  After  an 
acquire  cycle  the  data  is  dumped  on  to  a  standard  9 
track  computer  tape  as  a  single  block  headed  by  64  bytes 
of  Information  relating  to  current  channel  sample  rates, 
channel  memory  allocations  etc. 

The  12  analogue  Input  channels  are  converted  by 
Analogue  Devices  Corporation  A/D  converter  types  ADC10Z 
and  ADC10QU  with  20  and  8  vs  conversion  time  respectively. 

(These  convert  to  10  bit  accuracy  of  which  the  most  signi¬ 
ficant  8  were  used) .  Programing  the  system  to  set  channel 
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sample  rates  and  allocate  the  required  aeaory  areas  to 
the  channels  in  use  is  done  via  a  set  of  front  panel 
switches.  Toward  the  latter  part  of  the  experiment , 
data  not  specifically  requiring  the  attributes 
of  the  system  described  above  was  transferred  direct 
to  the  computer  memory  from  the  experiment. 

In  this  alternative  system,  based  on  the  Digital 
Equipment  Corporation  AR-11  multiplexing  A/D  converter 
module  with  crystal  controlled  real  time  clock  Inter¬ 
faced  with  the  Unibus  on  the  PDP-11/34  computer  (see 
Section  2.6),  the  maximum  throughput  rate  was  40kHZ, 
sufficient  for  most  purposes. 

Probe  transverslng  through  the  flow  was  accomplished 
with  a  stepper  motor  driven  linear  ball  screw  device 
moving  the  probe  carrier  In  0.0005" steps  at  variable 
rates  up  to  5.5kHZ.  Clock  pulses  from  the  traverse 
motor  control  were  used  to  drive  the  data  acquisition 
system  sample  conaands,  ensuring  perfect  synchronism 
between  the  probe  position  and  the  sequential  location 
of  the  measured  variable  in  the  output  data  record. 

2. 3  Optics  and  Photography 

Part  of  the  motivation  for  using  the  ozone-nitric 
oxide  reaction  was  the  separability  of  the  optical 
absorption  spectruras  (Figure  2.4).  This  property 
enables  either  ozone  or  nitrogen  dioxide  to  be  detected 
Independently,  but  it  was  possible  to  exploit  this  even 
further  by  taking  simultaneous  shadowgraph  pictures  of 
these  components  in  the  flow,  separating  the  two  wave- 


Fig.  2.4  Absorption  Spectra  of  0^.  anci  NO,. 
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length*  after  passing  through  the  test  section  with  a 
dlchrolc  beam  splitting  assembly. 

The  use  of  s  flla  piste  sised  Interference 
dlchrolc  beaa  splitter  and  Interference  filters  could 
be  avoided  by  using  Instead  a  combination  of  small 
narrow  band  pass  filters  at  the  light  source,  and  more 
economical  wide  band  components  at  the  laage  end.  The 
narrow  band  source  filters  were  custom  made  and  had  the 
following  characteristics. 

Centre  Wavelength  Peak  Transmission  50%  bandwidth 
O.V.  250  no  30%  ±  50  nm 

Blue  400  no  60%  t  50  nm 

At  the  Image  end,  the  bean  splitter  was  a  piece  of 
150  am  square  x  3  mb  thick  optical  quality  glass  plate 
coated  with  aluminium  to  75%  reflection  at  45°  in  the 
visible  region.  The  3  an  substrate  was  completely  opaque 
to  all  UV  coming  from  the  light  source.  Blue  light  was 
filtered  out  from  the  reflected  beam  by  a  150  x  75  mm 
optically  polished  conaercial  grade  "Short  Have  Filter" 
normally  used  In  light  sources  for  activating  fluorescent 
materials.  This  had  a  peak  transmission  wavelength  near 
250  nm  and  its  band  extended  to  the  near  ultraviolet 
allowing  some  contamination  of  the  UV  image  by  the  blue 
image.  However,  the  plate  required  over  8  exposures  by  the 
blue  source  to  achieve  the  same  exposure  as  a  single  exposure 
from  the  UV  source  which  is  quite  acceptable  for  the  current 
purpose.  The  layout  is  shown  in  2.5.  For  a  review  of 
ultraviolet  technology  the  reader  is  referred  to  Green  (1966). 
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An  air-arc  pulsed  light  source  emitted  e  wide 
spectral  band  extending  across  the  range  250  -  400  «■* 
being  essentially  a  high  temperature  black  body 
radiator.  A  5.5  uP  capacitor  provided  about  44J  of 
energy  to  the  spark  gap  giving  a  pulse  with  a  half  power 
duration  of  about  1.5  microseconds.  Ultraviolet  grade 
silica  glass  windows  enclosed  the  test  section  and 
aluminium  coatings  gave  adequate  reflectivities  at  UV 
and  blue  wavelengths.  Great  care  was  taken  to  prevent 
nitrous  oxides  contacting  the  reflective  surfaces  at 
any  time.  A  thin  coating  of  magnesium  flourlde  was 
applied  to  the  beam  splitter  to  protect  It  from  accidental 
exposure  and  long  term  degradation.  Standard  5"  x  7" 
flat  film  negatives  (Ilford  HP4)  were  found  to  have 
sufficient  UV  response  when  developed  to  about  1000ASA 
in  Ilford  Microphen  developer. 

2.4  Ozone  System 

The  unique  properties  of  ozone  demand  special  care 
and  handling  techniques.  It  Is  a  highly  unstable  allotrope 
of  oxygen,  making  It  a  particularly  powerful  oxidising  agent 
which  easily  reverts  to  diatomic  oxygen  on  contact  with 
material  surfaces.  Ozone  concentration  decays  with  time 
at  a  rate  proportional  to  the  surface  area  to  volume  ratio 
and  the  material  of  its  container,  with  rubber,  fabrics 
and  most  plastics  causing  a  faster  decay  than  metal  or 
glass*  (Sobersky  and  Slmena,  1973). 

t  The  only  plastics  used  in  the  apparatus  were  polytetra 

fluorethylene  (PTFE) ,  fluorlnated  ethylene  polythene  (FEP), 
Perspex,  and  polyester  resin.  None  of  these  exhibited  un¬ 
due  deterioration  after  repeated  exposures  to  high  concentra¬ 
tion  ozone. 


Added  Co  its  passive  decay  behaviour  is  its  very  real 
propensity  co  violent,  spontaneous  detonation,  being 
"■uch  more  sensitive  toward  detonation  than  hydrogen-oxygen 
aixtures"  (Harper  and  Gordon,  1959).  Harper  and  Gordon 
found  that  the  lower  concentration  limit  for  detonation 
of  an  oxygen-ozone  mixture  was  only  9Z,  at  which  point  a 
shock  wave  with  a  temperature  of  400°C  was  adequate  for 
ignition.  It  is  also  commonly  known  (and  demonstrated 
within  the  present  rig)  that  unexpected  explosions  occur, 
presumably  Initiated  by  the  presence  of  small  quantities 
of  organic  Impurities  in  the  apparatus.  The  requirements 
for  cleanliness  become  more  stringent  as  the  ozone  con¬ 
centration  Increases.  Naturally  the  Incompatibility  of 
ozone  with  organic  tissues  makes  it  very  noxious  to  the 
human  respiratory  system  as  well:  it  is  well  documented 
as  an  atmospheric  pollutant.  Interfering  with  oxygen  access 
in  the  lungs.  All  of  these  factors  were  considered  when 
designing  the  ozone  system. 

The  main  requirement  of  the  ozone  system.  Figure  2.6, 
was  to  as  safely  and  conveniently  as  possible  provide  up  to 
around  10%  ozone  in  a  passive  carrier  gas.  Ozone  was 
produced  by  a  conventional  electric  discharge  ozonator 
(Wallace  and  Tieman,  Model  BA023)  at  concentrations  of 
1  to  2%  in  oxygen.  It  was  then  separated  from  the  oxygen 
by  a  low  temperature  sorption  process  described  by  Cook 
at  al  (1966).  Approximately  1  kg  of  chromatographic  grade 
(28-200  mesh)  silica  gel  was  contained  in  a  stainless  steel 
pressure  vessel  which  was  placed  inside  a  large  dewar  flask 
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Fig.  2.6  Ozone  Generation  and  Concentration  System 


above  a  reserve  of  liquid  nitrogen.  A  1.5  on  thick 
cylindrical  copper  sheath,  laid  against  the  Inside  of 
the  devar  and  dipping  into  the  liquid  nitrogen,  helped 
distribute  the  temperature  Inside  the  flask.  This 
sheath  was  experimentally  shaped  and  fenestrated  to 
maintain  a  temperature  In  the  silica  gel  somewhere  be¬ 
tween  -  150°C  and  -  100°C.  Overcooling  was  found  to 
cause  condensation  problems  in  the  silica  gel,  and  at 
above  -  80°C  the  adsorption  capacity  of  the  gel  was 
insufficient.  The  dry  ozone-oxygen  mixture  from  the 
ozonator  passed  through  this  cell,  which  absorbed  ozone 
in  preference  to  oxygen.  After  the  required  quantity 
of  ozone  had  been  loaded  into  the  silica  gel,  the  dewar 
flask  was  removed  and  a  water  jacket  placed  around  the 
cell  as  the  carrier  gas  (nitrogen,  argon,  etc)  was  drawn 
through  to  pick  up  the  liberated  ozone  and  carry  it  to 
one  of  the  FEP  storage  bags  in  the  blow  down  rig. 

The  reliability  of  the  process  is  assured  when 
several  experimental  maxims  are  followed.  Firstly,  only 
low  decay  constant  materials  must  contact  the  ozone,  and 
they  must  be  scrupulously  clean.  In  this  rig,  the  only 
materials  were  stainless  steel,  aluminium  and  glass  as 
well  as  the  plastics  mentioned  above,  and  all  surfaces 
were  cleaned  with  xylene,  acetone  and  methanol  (Shea,  1976) 
before  assembly.  Secondly,  it  was  beneficial  to  ensure 
a  good  porosity  in  the  silica  gel  before  each  adsorption 
by  temporarily  fluidising  it  with  a  burst  of  dry  oxygen. 

A  mechanical  stirrer  was  particularly  useful  while  desorbing 
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since  It  assisted  heat  transfer  and  gas  flow  within 
the  cell.  The  paddle  of  the  stirrer  was  helical  and 
rotated  In  the  sense  to  lift  the  silica  gel.  All  other 
forms  tried  were  found  to  compact  the  gel  and  jam, 
particularly  at  low  temperature.  Thirdly,  the  oxygen 
supply  must  be  kept  as  dry  as  possible  to  avoid  blockage 
at  the  cell  by  ice  crystals,  and  to  preserve  the  adsorbtlve 
capacity  of  the  silica  gel.  It  was  necessary  to  regenerate 
the  silica  gel  after  approximately  every  twenty  cycles  by 
heating  in  air  to  160°C  for  2  to  4  hours. 

The  risk  of  exposing  laboratory  personnel  to  ozone 
was  minimised  by  forced  ventilation  in  the  vicinity  of 
the  rig  and  passing  all  vents  from  the  apparatus  through 
an  ozone  catalyst.  Both  soda  lime  and  activated  charcoal 
were  effective  for  this  purpose.  The  whole  system  from 
the  check  valve  after  the  ozonator  to  the  delivery  valve, 

was  designed  to  withstand  at  least  14  MPa,  a  worst  case 

detonation  condition. 

2.5  Instrumentation 
2.5.1  Velocity  Measurement 

The  average  velocity  in  the  x  direction  was  obtained 
from  both  the  pitot  static  differential  and  the  local 
average  temperature  using  the  relation. 

u  -  (2RT  (Pg  -  PJ/P)** 

where  P  -  static  pressure 

Pg  ■  stagnation  pressure 

T  *  local  time  mean  temperature 


R  ■  gas  constant 
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For  the  local  temperature  measurement  technique,  see 
2.5.2.  A  Datametrics  type  1014A  fast  response  differential 
electronic  manometer  was  connected  to  the  pitot  static 
tube  assembly  detailed  in  Figure  2.7(a).  The  signal  was 
recorded  in  digital  form  and  subsequently  smoothed  with 
a  low  pass  zero  phase  shift  recursive  digital  filter 
described  in  Appendix  C. 

2.5.2  Temperature  Measurement 

In  the  earlier  stages  of  the  project,  a  tungsten 
wire  operated  in  the  constant  current  mode  (l.e.,  as  a 
resistance  thermometer)  was  traversed  across  the  shear 
layer  to  obtain  a  mean  temperature  profile.  The  probe 
detailed  in  Figure  2.7(a),  was  made  from  two  0.3  mm 
diameter  steel  needles  epoxled  into  a  length  of  5  micron 
diameter  hollow  ceramic  rod.  A  piece  of  3  micron  tungsten 
wire,  copper  plated  except  for  the  middle  1  mm,  was  soft 
soldered  across  the  prongs.  The  resulting  active  portion 
of  the  wire  therefore  had  a  length-to-diameter  ratio  of 
about  200,  and  a  cold  resistance  of  about  7  ohms.  (These 
wires  were  used  by  Thomas  (1977)  as  hot  wire  anemometers). 

A  simple  circuit  (Figure  2.7(b)  )  provided  a  1.0  mA 
constant  current  to  measure  the  wire  resistance.  An 
energy  balance  on  the  wire  in  still  air  shows  that  the 
temperature  rise  due  to  electrical  heating  is  less  than 
0.2°R,  which  is  negligible  compared  to  the  fluid  tempera¬ 
ture  variations  and  hence  the  probe  calibration  is 
velocity  independent.  Neglecting  the  0.2°K  overheat,  the 
small  higher  order  non-linearities  in  the  temperature 
coefficient  of  resistivity,  and  assuming  a  perfectly  con- 
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etant  current,  I,  it  Is  easily  shown  that 
AV/AT  -  1  «  Rf 

where  is  the  circuit  feedback  resistance  and  «  is 
the  temperature  coefficient  of  resistivity  of  the  wire 
material.  The  quantity  on  the  right  is  easily  maintained 
constant  within  the  range  of  experimental  error,  so  adequate 
linearity  is  assured. 

Since  the  calibration  was  somewhat  subject  to  drift, 

(see  below) ,  it  was  necessary  to  clean  and  calibrate 

it  prior  to  each  tun.  Saturated  steam  was  used  as  a  reference 

temperature. 

The  frequency  response  is  that  of  a  first  order 
linear  system,  with  a  cutoff  frequency  around  100  Hz  in 
air  at  15  m/s. 

This  type  of  probe  had  several  drawbacks  in  use,  which 
became  more  severe  with  higher  temperature  runs.  The  nit¬ 
rous  chemicals  present  in  the  test  section  were  highly 
hygroscopic,  leading  to  deposition  of  acidic  droplets  on 
the  probe,  which  in  turn  led  to  rapid  corrosion  and  unstable 
resistance.  Besides  this,  the  soldered  junctions  at  each 
end  of  the  wire  produced  some  component  of  thermovoltalc 
asymmetry  which  varied  from  probe  to  probe. + 

In  order  to  make  reliable  temperature  measurements  in 
rims  exceeding  about  30°K  rise,  it  was  necessary  to  develop 
a  fast  thermocouple  probe  (Figure  2.8(a)  ).  A  chromel- 


t  The  simplest  way  to  check  was  to  move  the  probe  through  a  jet 
of  hot  air  and  watch  the  recovery  of  the  amplifier  output. 

The  thermovoltalc  component  has  a  comparitlvely  long  time 
constant. 
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Fig.  2.8(a)  Detail  of  Thermocouple  Probe. 


Fig.  2.8(b)  Thermocouple  Amplifier. 


alumel  wire  pair,  13  micron  diameter,  was  welded  together 
to  form  the  active  junction.  These  were  soft  soldered 
to  0.25  mm  prongs  of  the  same  material  anchored  In  a 
piece  of  3  mm  diameter  dual  bore  ceramic  Insulator.  The 
frequency  cutoff  for  this  probe  was  measured  at  around 
200  Hz  in  air  at  20  m/s  which  was  quite  adequate  to 
resolve  a  mean  temperature  profile  (Figure  2.9).  As 
with  the  resistance  wire,  the  sensor  was  prone  to  hygro¬ 
scopic  formation  of  droplets,  but  these  would  not  change 
the  calibration  unless  they  bridged  the  actual  welded 
junction,  which  could  easily  be  prevented  by  washing  in 
distilled  water  and  drying  immediately  before  a  run. 
Corrosion  was  still  a  problem,  which  could  be  minimized 
by  storing  over  dry  silica  gel,  but  no  evidence  of  slgnifl- 
cant  calibration  error  could  be  found  to  the  point  where 
the  probe  went  open  circuit.  The  mean  life  was  4  or  5 
runs.  Internal  resistance  of  the  thermocouple  with  leads 
was  typically  45  ft,  and  this  was  taken  into  account  when 
calculating  the  voltage  gain  of  the  circuit  shown.  (Fig¬ 
ure  2.8(b)  ).  Since  a  temperature  rise  was  being  measured, 
there  was  no  necessity  for  a  reference  junction. 

2.5.3  Pressure  Gradient 

The  static  pressure  gradient  was  measured  in  the 
low  speed  stream  with  two  static  pressure  probes  (similar 
to  the  one  in  Figure  2.6)  placed  at  X  ■  0  and  100  mm.  The 
pressure  was  measured  by  a  Da tame  trie  electronic  manometer 
and  differential  pressure  head  and  the  signal  filtered  with 
an  RC  filter  with  a  50  Hz  cutoff  frequency.  Initial  runs 
proved  the  pressure  gradient  to  remain  constant  throughout 
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Fig.  2.9  Raw  Temperature  Probe  Data  from  Traverse  of  Reacting  Shear 
Layer. 
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Che  run,  Chun  requiring  only  a  low  sample  rata,  and 
this  filtering  removed  the  rlak  of  error  due  to  Nyqulst 
folding  about  the  sample  frequency. 

2.5.4  Ozone  Concentration 

Due  to  the  unstable  nature  of  ozone  and  the  number 
of  variables  in  the  ozone  generating  process,  it  was 
necessary  to  measure  the  ozone  concentration  immediately 
upstream  of  the  test  section  as  indicated  in  Figure  2.2. 
The  ultavlolet  absorption  meter  is  detailed  in  Figure 
2.10.  It  operates  on  the  principle  of  Beer's  Law. 

viz  I/I  -  e_kn°3d 

o 

where  Iq  ■  incident  light  intensity 

I  ■  transmitted  light  intensity 

k  -  absorptivity 

no  3  ■  ozone  number  density 
d  ■  path  length 

The  construction  of  the  device  made  it  possible  to  set 
the  gap  with  micrometer  accuracy.  The  light  source 
was  a  Penray  Model  SC-1  low  pressure  mercury  lamp  which 
is  very  stable  and  radiates  most  energy  at  253.7  mm.  It 
was  driven  by  a  regulated  inverter  power  supply  running 
at  20  kHz  (to  make  it  inaudible).  Two  optical  filters 
were  placed  in  series  to  isolate  the  253.7  mm  mercury 
line.  One  was  an  absorption  filter  supplied  by  the 
Penray  Company  and  the  second  was  a  bandpass  interference 
filter  with  a  50Z  intensity  bandpass  of  240  -  260  nm. 
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Fig.  2.10(a)  Arrangement  of  Ozone  Monitor 
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Having  effectively  Isolated  the  single  wavelength  . 
the  value  of  k  was  known  to  extremely  good  eccuracy 
(Griggs,  1968). 

A  photodetector  supplied  by  United  Detector 
Technology,  model  Photop  500  with  special  UV  en¬ 
hancement  was  used  to  measure  the  transmitted  light 
intensity.  The  circuit  detailed  is  a  simple  trans¬ 
conductance  amplifier  with  provision  for  offsetting 
the  output  to  zero  volts  under  dark  conditions,  and 
an  RC  filter  to  remove  noise  and  ripple  from  the  20  kHz 
lamp  supply. 

The  lamp  was  turned  on  at  least  30  minutes  before 
a  run  to  give  it  time  to  stabilize.  Immediately  before 
the  run,  the  lamp  was  shielded  and  the  output  adjusted 
to  zero  with  the  offset  control.  The  output  was  sampled 
and  digitized  immediately  before  the  solenoid  operated 
valves  were  given  the  open  command  in  order  to  establish 

the  value  of  V  . 

0 

Since  the  system  did  not  use  a  reference  beam,  the 
drift  rate  was  very  carefully  checked,  and  found  to  be 
less  than  0.7%  in  V0  and  less  than  0.11  in  Vj  over  a 
period  of  5  minutes. 

2.6  Computing  Facilities 

Throughout  the  study  a  PDP11/34  processor  was  used, 
-with  28k  of  volatile  memory.  Peripheral  devices  were  two 
RK05  disk  drives,  an  IBM  nine  track  tape  drive,  a  line 
printer,  a  teletype  terminal  and  a  Hewlett  Packard  7221A 
plotter. 


3.  RESULTS 


3.1  Photographic  Investigation 

Due  to  the  unusual  nature  of  this  flow  visualization, 
this  study  will  be  more  meaningful  if  the  optical  properties 
of  the  system  are  recognized.  The  information  seen  on  the 
photographs  is  in  general  some  combination  of  refractive 
index  information  and  species  absorption  Information.  The 
refractive  index  information  is  the  conventional  shadow¬ 
graph,  where  the  angular  deflection  of  a  ray  is  a  function 
of  the  normal  refractive  index  gradient  integrated  along 
its  path.  Refractive  index  is  a  function  of  denalty,  wave¬ 
length  and  species,  so  both  thermal  gradients  and  species 
gradients  (when  present)  will  appear  on  the  plates.  The 
refractive  index  signature  must  similarly  appear  in  both 
UV  and  blue  shadowgraphs,  but  due  to  dissimilar  light  in¬ 
tensity  as  well  as  film  plate  characteristics  the  contrast 
is  different.  Due  to  the  large  optical  path  length  between 
the  test  section  and  film  plate  (Figure  2.5)  the  bent  ray 
becomes  considerably  displaced  by  the  time  it  reaches  the 
film,  which  accounts  for  some  loss  in  resolution.  The 
light  absorption  by  a  gaseous  species  is  an  exponential 
function  of  density,  path  length  and  absorptivity.  The 
optical  density-exposure  characteristic  of  the  film  is 
logarithmic  so  that  the  optical  density  at  a  point,  within 
a  certain  range,  is  a  linear  function  of  the  total  absorb¬ 
ing  species  within  the  ray.  At  the  wavelengths  used,  how¬ 
ever  there  is  e  factor  of  20  difference  between  the  absorp- 
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tion  coefficients  of  nitrogen  dioxide  and  ozone.  The 
mean  concentration  required  to  register  a  difference 
of  one  stop  on  the  plate  is  O.SZ  and  0.02Z  respectively, 
so  the  sensitivity  of  the  UV  plate  to  ozone  is  far 
greater  than  the  blue  plate  to  nitrogen  dioxide.  Effects 
from  the  boundary  layer  on  the  viewing  windows  are  general¬ 
ly  negligible  except  for  the  ozone  pictures.  Ho  signifi¬ 
cant  refractive  bending  would  occur  because  the  boundary 
layers  are  thin  (-  0.5  on)  and  the  gradients  tend  to  be 
aligned  with  the  light  beam.  However  in  the  case  where 
the  boundary  layers  are  locally  contaminated  with  a 
typical  maximum  possible  concentration  of  nitrogen  dioxide 
or  ozone,  the  appropriate  light  beam  would  be  attenuated 
to  the  extent  of  0.1  and  2  stops  respectively.  Some  tests 
were  done  to  determine  how  good  the  dlchrolc  optics  were 
in  separating  the  frequencies.  Contamination  of  the  UV 
plate  by  blue  light  was  found  to  be  3  stops  down,  while 
UV  contamination  of  che  blue  plate  was  immeasurable. 

Smear  due  to  finite  exposure  time  is  insignificant, 
amounting  to  about  50  urn  (.002")  on  the  full  size  plate. 
Finally,  all  pictures  were  taken  through  a  100  nan  width 
of  flow,  about  4  of  the  maximum  visual  shear  layer  thick¬ 
nesses  seen  in  the  photographs.  In  discussiog  the  photo¬ 
graphs,  for  brevity  the  location  along  the  layer  will  be 
stated  simply  as  the  proportion  of  the  total  length  from 
the  left  hand  edge,  e.g.,  x  ■  l  refers  to  the  extreme 
right  hand  edge. 

The  first  set  of  pictures  is  from  a  series  of  runs 
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In  the  He/N^  shear  layer.  In  which  the  heat  release  was 

progressively  Increased  while  keeping  all  other  variables 

constant  (Figures  3.1  -  3.4).  Figure  3.2  is  from  an 

identical  flow  as  3.1.  It  is  included  to  demonstrate  the 

variability  between  different  realisations.  The  Reynolds 

number  at  x  -  1  in  the  non-reacting  pictures  (3.1,  3.2)  is 
Au6 

— -  4  x  104  and  since  the  maximum  temperature  that  can 

N2 

exist  at  the  flame  sheet  is  300  and  420°C  in  Figures  3.3 
and  3.4  respectively,  the  corresponding  local  Reynolds 
numbers  in  these  regions  of  maximum  temperature  are  de¬ 
creased  to  352  and  27Z  of  the  cold  value.  A  detailed 
examination  of  these  three  pictures  reveals  several  effects 
as  a  result  of  heat  release.  Firstly,  the  initial  roll-up 
of  the  Kelvln-Helmholtz  instability  in  the  laminar  layer 
appears  to  be  delayed.  Secondly,  there  is  a  dramatic 
change  in  the  appearance  of  the  small  scales.  Heat  release 
apparently  has  the  effect  of  delaying  transition  and  in¬ 
hibiting  the  energy  cascade  to  higher  frequencies.  It  is 
worth  noting  here  that  the  Kolmogorov  microscale  varies 
as  T1’3  but  is  well  below  the  resolution  of  the  photographs. 

The  behaviour  of  the  small  scales  would  appear  to  be  some 
function  of  the  local  Reynolds  number  (based  on  the  temp¬ 
erature  in  the  local  flame  sheet).  Clearly  the  dynamics 
of  these  small  scales  must  be  directly  effected  by  volume 
expansion  as  well  as  simple  viscosity  increases,  and  this 
may  be  the  dominant  factor  in  the  suppression  of  the 
highest  frequencies.  Hence,  the  effect  of  heat  release 
in  the  shear  layer  is  a  stabilising  one,  since  the  moleculer 
scsle  mixing  processes  upon  which  the  resctlon  itself  ultimately 


FIGURE  3.2  Shadowgraph  of  Non-Reacting  He/N  Shear  Layer 
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FIGURE  3.3 
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depends  appear  to  become  less  efficient.  (Obviously 
it  would  be  questionable  to  draw  such  conclusions 
merely  from  the  three  realizations  shown;  however, 
this  observation  is  generally  consistent  with  a 
number  of  other  realizations,  including  those  made 
in  other  gases.  The  main  reason  for  exhibiting  the 
heliian-nltrogen  layer  here  is  that  much  more  information 
is  present  due  to  the  refractive  index  difference). 

A  priori  it  is  not  clear  whether  heat  release 
increases  or  decreases  the  width  of  the  layer.  Certain* 
ly  the  volume  expansion  due  to  heat  release  causes  a 
v  component  of  velocity  at  infinity.  However  there  is 
no  significant  width  increase  indicated  on  the  pictures. 
Measurements  presented  below  will  confirm  this.  Since 
the  density  is  lower,  the  amount  of  mass  within  the 
layer  is  therefore  decreased.  Implying  a  reduced  overall 
entrainment  of  free  stream  fluid.  The  extent  of  en¬ 
trainment  reduction  is  apparently  enough  to  more  than 
offset  the  mean  volume  expansion  within  the  layer.  From 
the  photographs  it  is  not  clear  what  effect  the  heat 
release  has  had  on  the  dynamics  of  the  large  structure 
that  might  account  for  this  greatly  reduced  entrainment. 
Nevertheless  the  presence  of  a  large  structure  remains 
clear  at  least  up  to  the  temperature  rises  and  Reynolds 
numbers  in  this  experiment. 

A  series  of  photographs  of  nitrogen-nitrogen  shear 
layers  was  also  taken  (Figures  3.5  to  3.8),  using  the 
same  velocities  as  for  the  previous  set.  The  maximum 
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mean  temperature  rises  are  4.4,  77,  98,  and  122°C 

respectively.  Since  the  refractive  indexes  of  .the 
free  stream  gases  are  the  sane  only  the  thermal 
gradients  (not  the  species  gradients)  will  be  visible. 

For  this  reason,  and  due  to  the  high  level  of  small 
scale  turbulence,  large  scales  in  the  layer  are  not  as 
clear  as  in  the  helium-nitrogen  layer.  However,  there 
is  evidence  of  large  scale  clumping  of  the  nitrogen 
dioxide  in  Figure  3.8.  These  nitrogen-nitrogen  re¬ 
acting  shear  layers  show  similar  effects  of  heat 
release  as  the  nitrogen-helium  shear  layers.  The 
overall  growth  rate  changes  very  little,  and  the 
finest  turbulence  scales  are  progressively  suppressed. 

looking  at  the  UV  pictures  of  ozone  associated  with 
all  the  above  figures,  it  is  clear  that  the  penetration 
of  ozone  across  the  layer  depends  on  the  reactant  con¬ 
centration  ratio  £,  and  the  level  of  small  scale  turb¬ 
ulence.  (It  must  be  remembered  that  the  frequency  of 
UV  light  used  was  highly  sensitive  to  ozone,  and  that 
any  mean  ozone  concentration  above  about  0.2%  will  be 
opaque).  The  engulfment  process  is  particularly  clear 
in  Figure  3.3  at  x  ■  0.5.  The  ozone  carrying  stream  is 
entwined  into  the  periphery  of  the  clump  of  vorticity, 
and  carried  around  by  it  while  reaction  is  occuring  at 
the  interface  between  the  two  gas  streams.  This  appears 
to  be  occuring  mainly  by  laminar  diffusion  from  both  sides, 
and  is  seen  to  thin  down  the  "tongue"  of  ozone  until  it 
disappears  entirely  before  achieving  %  a  rotation  of  the 
large  structure.  The  effect  of  introducing  some  small 
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scale  turbulence  can  be  observed  by  contrasting 
this  with  a  structure  further  downstream  at  x  -  0.8, 
where  the  "tongue"  of  ozone  is  much  more  rapidly  con¬ 
sumed,  particularly  from  the  inside  edge  where  there  is 
a  higher  level  of  turbulence.  A  similar  behaviour 
is  seen  in  the  nitrogen-nitrogen  shear  layer  pictures 
(e.g. ,  Figure  3.6)  but  the  small  scale  turbulence  is 
even  more  energetic  and  results  in  a  comparatively 
short  life  of  entrained  ozone.  Clearly,  Increasing 
the  ozone  concentration  relative  to  nitric  oxide 
concentration  causes  ozone  to  penetrate  further  across 
the  layer,  and  in  the  case  of  Figure  3.8,  ozone  has 
become  the  majority  species  (l.e,  there  is  an  excess  of 
ozone  after  all  the  nitric  oxide  Is  burnt)  almost 
throughout  the  layer.  These  pictures  provide  a  qualita¬ 
tive  support  for  the  PDFs  measured  by  Konrad  (1977), 

(see  Appendix  Al),  if  it  is  borne  in  mind  that  mixtures 
of  mixture  fraction  c  where  0  £  c  <  cgappear  transparent 
to  the  UV  light,  and  mixtures  where  cg  <  c  ^  l  appear 
opaque. 

The  mechanics  of  the  large  structure  under  conditions 
of  heat  release  can  be  usefully  studied  prior  to  transi¬ 
tion.  The  growth  rate  and  entrainment  rate  are  evidently 
quite  similar  before  and  after  transition,  the  predominant 
change  being  an  Increase  in  the  molecular  scale  mixedness 
(Konrad  1976,  Breidenthal  1978).  Figure  3.9  shows  two 
flows  with  temperature  rises  of  10  and  160°C  respectively. 
The  maximum  Reynolds  number  in  each  case  is  about  10u. 


FIGURE  3.9  Blue  Shadowgraphs  of  Low  Speed  Reacting  N  /N 

Shear  Layers 
Upper  Stream  5  ms'1 
Lower  Stream  1  ms'1 
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The  absence  of  fine  scale  turbulence  permits  a  clear 
view  of  a  dramatic  change  In  the  character  of  the 
large  structures.  A  measure  of  the  relative  strength 
of  vortielty  clumping  can  be  seen  In  the  width  of  the 
connecting  braid  between  adjacent  structures,  and  by 
this  measure  alone  It  Is  observed  that  the  vortielty 
clumping  Is  weakened  with  Increasing  reaction.  Al¬ 
though  the  far  field  strength  of  the  induced  velocity 
field  around  an  isolated  patch  of  vortielty  depends  on 
the  total  vortielty,  not  on  how  it  Is  distributed,  the 
induced  velocity  in  the  near  field,  being  a  summation 
effect  of  distributed  vortielty  will  be  somewhat  less. 

This  will  tend  to  inhibit  two  of  the  mechanisms  res¬ 
ponsible  for  entraining  free  steam  fluid  into  the  layer. 

The  first  is  the  entwining  process,  causing  a  particular 
eddy  to  grow  by  simple  acquisition  of  fluid  and  the  second 
is  the  vortex  pairing  process,  which  is  observed  to  en¬ 
train  by  engulfment  of  free  steam  fluid. 

A  further  effect  on  entTainment,  present  after  tran¬ 
sition,  is  that  due  to  the  wiping  out  of  the  higher 
frequency  turbulence.  However  this  is  a  small  component 
of  the  total  entrainment  (Konrad  reports  eight  percent). 
3.1.1  The  Role  of  Baroclinlc  Torque 

The  possible  action  of  baroclinlc  torque*  in  causing 
non-reacting  shear  layers  with  unequal  free  stream  densities 
to  have  unequal  entrainment  rates  from  each  of  the  free 


t  Many  of  the  ideas  in  this  section  arose  during  discussions 
with  Professor  G.L.  Brown. 
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streams  is  discussed  In  section  3.5.  The  barocllnlc 
torque  action  Is  used  to  show  how  the  vorticity  Is 
transferred  to  the  less  dense  side  of  the  layer  and 
thereby  Induces  a  greater  proportion  of  less  dense 
fluid  to  enter  the  layer.  Barocllnlc  torque  is  a  re¬ 
sult  of  non-aligned  density  gradient  and  acceleration 
vectors,  which  are  obviously  a  feature  of  mixing  flows 
with  chemical  heat  release.  The  derivation  of  the 
expression  for  barocllnlc  torque  from  Euler's  equation 
Is  relatively  simple,  but  summarized  here  for  the  sake 
of  completeness. 

Euler's  equation 

3u  , 

^  +  (a  •  7)  u  *  -  -  ?  p  +  g 

Taking  the  curl,  and  neglecting  gravity 

3<d  Vp 

j”  +  V  x  (u*  V)  u  ■  -  (V  x  ~) 


3u  . 

■“■  +  (u  •  V)  id  ■  (V  p  x  Vp) 


substituting  for  Vp  from  Euler's  equation  gives 


Did 

Dt 


1  °U 
p  (!P  X  Dt  ) 


That  barocllnlc  torque  is  relevant  to  the  question 
of  entrainment  reduction  in  the  reacting  shear  layer 
depends  on  the  view  that  the  mechanics  of  the  layer  are 
controlled  by  large  structure  which  is  basically  two 
dimensional.  Some  idea  of  the  distribution  and  action 
of  barocllnlc  torque  can  be  obtained  by  considering  the 
mean  density  gradient  and  acceleration  vector  around  the 
large  structure  shown  In  Figure  3.9.  It  Is  shown  In 
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Appendix  A  that  MO  concentration,  as  seen  in  Figure  3.9, 

2 

is  analogous  to  temperature  rise  (or  density  decrease) 

since  the  Schmidt  and  Lewis  numbers  are  approximately 

one.  An  envelope  drawn  along  the  Interface  of  this  hot 

(turbulent)  region  with  the  non-turbulent  free  stream 

fluid  then  defines  a  surface  normal  to  the  mean  density 

gradient.  Such  a  surface  is  shown  schematically  as  a 

dashed  line  in  Figure  3.10,  superimposed  on  a  representation 

of  the  mean  streamlines  (where  the  observer  is  moving 

with  the  large  structure).  The  flow  pattern  shown  is 

similar  to  that  given  by  Patnaik,  Sherman  and  Corcos 

(1976).  Careful  inspection  of  the  vector  product  of  Vo 
Du 

and  around  the  dashed  envelope  indicates  that  the 
dominant  contribution  to  the  vorticty  by  barocllnlc 
torque  action  occurs  at  the  points  labelled  'A' ,  and 
this  contribution  is  in  the  same  direction  as  the  general 
circulation.  Control  volume  analysis  shows  that  in  any 
given  length  of  the  layer,  the  mean  total  circulation 
must  remain  constant,  implying  that  the  vortlcity  at 
points  A  is  growing  at  the  expense  of  the  vorticity 
concentrated  in  the  core  of  the  large  structures  and/or 
in  the  braids.  This  action  is  qualitatively  consistent 
with  the  effect  of  heat  release  shown  in  Figure  3.9,  and 
provides  an  explanation  for  the  observed  reduction  in 
entrainment  rate  of  free  stream  fluid.  The  crucial  feature 
is  the  weakened  near  field  induced  velocity,  due  to  the 
more  diffuse  vortlcity.  In  particular,  components  of 
velocity  in  the  y  direction  will  be  relatively  weaker, 
thereby  slowing  both  the  entwining  and  engulfment  mechanisms 


streamlines 

edge  of 
turbulent  zone 


Fig.  3.10  Illustration  of  Baroclinic  Torque  Action 
Reacting  Shear  Layer  with  Heat  Release. 


slope  thickness.  Figure  3.11  shows  both  velocity  profiles 
and  y  ■  0  represents  the  dividing  streamline  calculated 
as  described  in  Appendix  B.  The  reasons  for  the  very 
slight  differences  in  profile  shape  are  not  clear  but  it 
seems  unlikely  to  be  due  to  some  effect  of  heat  release 
as  the  maximum  temperature  rise  was  only  20°C.  They 
have  been  plotted  together  to  confirm  the  similarity 


form  of  the  profiles  and  particularly  to  clarify  the 

relative  positions  of  the  temperature  profiles  shown  on 

the  same  figure.  The  typical  result  shown  exhibits  a 

widened  and  rightward  shifted  inferred  temperature  profile. 

Inspection  of  a  series  of  runs  made  with  varying  £  and 

AT  „  confirmed  this  trend.  In  a  sample  of  seven  of  these 
max 

results,  the  mean  shift  of  the  position  of  AT  was  0.26  , 

r  max  w 

and  the  standard  deviation  in  the  sample  was  0.116^. 

For  the  purposes  of  comparing  widths  of  mean  temperature 


profiles,  if  6^,  a  half  maximum  amplitude  width  is  defined, 
then  the  mean  value  of  6^/6^  in  the  above  sample  is  0.84, 
compared  with  1.26  for  the  inferred  temperature  profiles. 


The  value  of  AT 


is  used  later  as  a  parameter  for 


max 
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(x  -X<>) 


Fig.  3.11  Comparison  of  Velocity  and  Temperature  Profiles. 


for  comparison  between  experimental  and  Inferred 

results,  since  It  Is  a  readily  identifiable  point 

on  any  temperature  profile  which  is  directly 

dependent  on  the  entrainment  and  mixing  processes 

in  the  shear  layer.  Detailed  results  relating  to 

different  AT  in  the  experimental  and  inferred 
max 

profiles  are  presented  in  Section  3.4.2. 

It  is  not  clear  why  these  differences  between 
the  inferred  and  measured  temperature  profiles  exist 
It  is  possible  that  there  could  be  some  effect  of 
finite  reaction  rate  being  more  pronounced  on  the 
high  speed  side,  but  Increasing  the  reactant  con¬ 
centration  by  a  factor  of  five  (thereby  reducing 
Tchem  by  tbe  8ame  or<^er)  does  not  appear  to  change 
the  observed  behaviour.  There  is  some  uncertainty 
about  the  signal  to  noise  ratio  of  the  composition 
probe  used  by  Konrad  to  determine  the  PDFs  in  his 
uniform  density  flow+  consisting  of  a  mixture  of 
helium  and  argon  flowed  with  nitrogen  on  the  other 
side  of  the  layer.  This  was  not  a  problem  in  his 
helium-nitrogen  shear  layer  data  since  the  signal 
was  much  stronger.  However,  when  inferred  and 
measured  temperature  profiles  for  a  helium-nitrogen 
layer  are  plotted  (Figure  3.12)  the  discrepancy  is 
still  present.  For  this  figure  the  experimental 
profile  was  obtained  for  the  same  flow  parameters 
as  Konrad's  data  (r  -  0.38,  s  "  7.0)  and  since  the 

t  Private  communication.  Professor  G.L.  Brown. 
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Y/(X-XO) 


Fig.  3.12  Comoarison  of  Temnerature  Profiles. 
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the  dividing  streamline  could  not  be  determined 
from  this  experimental  data.  Its  position  was 
established  by  matching  the  pu7  profile  with 
Konrad's. 

3.2.2  Calculation  of  the  Mixture  Fraction  PDFs 

from  the  Experimentally  Measured  Temperature 
Profiles  -  the  Inverse  Problem. 

In  Appendix  A  it  Is  shown  how  the  temperature 
at  a  point  can  be  calculated  from  the  mixture  fraction 
PDF  and  the  free  stream  reactant  concentrations, 
given  the  heat  release  of  the  reaction,  the  specific 
heats,  and  the  assumption  of  infinite  reaction  speed. 
In  this  way  it  was  possible  to  infer  temperature 
rises  from  the  non-reacting  PDF  measured  by  Konrad, 
although  the  procedure  Is  strictly  only  valid  for 
negligible  temperature  rise  because  the  PDFs  could 
be  expected  to  change  under  heat  release  conditions. 
The  solution  for  the  mean  temperature  at  a  given 
location  is  simply 


■I. 


AT  «■  C0  N0>  ’  PC<C>T<C’  «•  C, 


NO 


)dc 


(3.1) 


where  £  ■  free  stream  reactant  concentration  ratio, 

C  /C  ,  and  p(c)  is  the  mixture  fraction  PDF. 

0  03  0  NO 

The  function  T  is  shown  in  Appendix  A  to  be 


AH 

t(c,  e,  c  )  -  c 

0  NO  0  NO  -p 


nf  r 

-~lc  H(cg  -  c)  +  £(1  -  c)H(c  -  cg) 

D  » 


where  cg  ■  C/(l  +  O,  the  stoichiometric  mixture  fraction, 
and  H(c)  is  the  Heaviside  step  function. 


Define 


T<c*  «•  co  NO)Cn 

T*  - - C - AH  “  cH(ce  “  c)  +  5(1  "  c)H(c  "  S) 

0  NO  f 

(3.2) 

(i.e. ,  scale  the  temperature  by  the  temperature  rise 
that  would  be  obtained  if  all  the  free  stream  HO  were 
burned) . 

Then 

AT*  (O  -  |pc(c)  T*  (C,  c)dc  (3.3) 

J0 

This  result  is  an  integral  equation  of  the  first  kind 
for  p(c)  if  AT*  (£)  and  T*  (£,  c)  are  known.  Although 
T*  is  known  everywhere  (equation  3.2),  each  run  of  the 
experiment  determines  AT*  for  only  a  single  value  of 
the  free  stream  reactant  concentration  ratio  €. 

It  is  clear  from  3.3  however  that  the  inert  PDF 
p(c)  could  be  obtained  experimentally  if  AT*  (O  is 
measured  for  a  wide  range  of  values  of  £.  The  equation 
could  be  inverted  in  a  number  of  ways,  one  of  which 
would  be  to  approximate  the  measured  AT*  (€)  by  a 
function  which  closely  fits  the  data  over  a  vide  range 
of  £.  This  method  of  determining  the  PDF  has  great 
merit  in  that  it  does  not  require  probes  with  the 
time  response  and  spacial  resolution  to  resolve  the 
Kolmogorov  microscale:  the  probe  is  only  required  to 
measure  the  time  mean  temperature  rise.  The  main 
difficulty  in  the  present  application  is  the  extent 
to  which  the  approximation  of  fast  chemistry  is  valid 
(see  Appendix  D). 
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The  function  p(c)  was  approximated  by  an  arbitrary 
mssber  n  delta  functions,  each  of  area  a^,  and  h  -  1 
values  of  AT*  at  various  £  were  chosen. 

n 

Pc<0  -  l  (3.4) 

Equation  3.3  was  then  solved  numerically  by  a  matrix 
method  to  solve  for  the  areas  of  n  -  1  of  the  delta 
functions,  the  last  being  found  from  the  PDF  property 
relation 

[  pc(c)dc  -  1 
j0 


lit*  ) 

n 

I  -  1  (3.5) 

i-1  1 

(The  nt(>  area  is  actually  the  sum  of  the  areas  of 
the  delta  functions  at  c  ■  0  and  c  ■  1:  this  inverse 
problem  cannot  solve  for  them  separately  because  at  c  »  0 
and  1,  T*  ■  0,  giving  these  areas  zero  coefficients  in 
the  matrix  equation).  This  attempt  to  solve  for  p(c) 
failed  because  the  quantity  and  precision  of  the  data 
were  insufficient  to  establish  the  correct  function  for 
AT*.  The  failure  criterion  was  that  pfi(c)  should  be 
positive  everywhere  in  the  range  0  <  c  <  1.  A  study 
of  the  disposition  of  the  matrix  equation  revealed  that 
it  would  be  more  sensitive  to  inaccuracies  in  AT*  as  n 
increased.  Consequently  a  value  of  n  *  5  was  chosen, 
because  of  all  the  experimental  data,  only  4  runs  fell 


within  the  acceptance  conditions  of  (a)  identical 
flow  velocities  and  gases  (b)  low  temperature  rise 
(  <  70°C)  (c)  wide,  well  spaced  range  of  £.  The 

necessity  of  fitting  a  function  AT*  was  thereby 
avoided.  The  procedure  was  repeated  for  AT*  taken 
at  seven  values  of  n,  thereby  creating  seven  PDFs 
for  the  mixture  fraction  at  these  locations  across 
the  layer.  Figure  3.13.  In  this  figure  it  must  be 
remembered  that  the  absolute  values  of  the  delta 
functions  corresponding  to  c  ■  0  and  1  are  not  de¬ 
termined,  only  their  sum  as  explained  above.  This 
area  has  been  arbitrarily  distributed  between  the  two 
pure  mixtures  in  the  same  ratio  as  the  mean  mixture 
fraction  resulting  from  a  combination  of  the  mixtures 
represented  by  the  other  four  delta  functions.  With 
the  exception  of  the  end  delta  functions.  Figure  3.13 
can  be  compared  with  Konrad's  PDFs  shown  on  the  next 
figure.  Although  the  details  in  the  PDFs  Inferred  by 
the  inverse  method  are  lacking,  it  is  apparent  that 
the  qualitative  and  quantitative  features  of  both  sets 
of  PDFs  are  compatible,  in  particular  the  evidence  of 
the  same  most  probable  mixture  being  present  all  across 
the  mixing  layer,  while  the  mean  mixture  fraction  is 
different  at  each  location. 

It  was  thought  that  by  assuming  a  hlstrogram  form 
for  Pc(c)  Instead  of  a  delta  function  form  would  lead 
to  a  physically  more  satisfactory  PDF.  This  form  also 
produced  regions  of  negative  probability  due  to  the 
measured  AT*  being  insufficiently  precise.  An  examination 


of  the  matrix  solution  procedure  indicated  that  a 
histogram  form  of  p£(c)  produced  relatively  email 
pivoting  coefficients  in  the  matrix,  making  the 
solution  very  sensitive  to  Inaccuracies  in  the  data. 
Further  refinements  did  not  seem  warranted  until 
data  for  a  more  comprehensive  set  of  values  of  £  are 
obtained. 

From  the  PDF  it  is  possible  to  estimate  c^,  the 
mean  mixture  fraction  in  the  turbulent  zone.  This  is 
obtained  by  taking  the  mean  mixture  fraction  of  the 
PDF  without  the  contribution  of  the  delta  functions  at 
c  ■  0  and  1.  A  plot  of  Cp  against  the  maximum  slope 
thickness  is  shown  in  Figure  3. 14  (a)  for  the  PDFs  of 
Konrad  and  the  present  work.  The  reason  for  the  steeper 
slope  of  the  present  result  is  not  clear,  but  it  is 
consistent  with  the  different  widths  of  the  temperature 
profiles  given  by  the  two  sets  of  PDFs  (see  Figure  3.11). 
In  fact,  when  the  plot  of  Cj  is  scaled  by  6^  instead  of 
6^,  the  agreement  is  remarkably  good  (Figure  3.14(b)  ) 
considering  the  vastly  different  methods  of  estimating 

c,p.  Despite  this  good  agreement,  the  problem  still 
remains  as  to  why  the  Cj  profiles  are  so  different 
relative  to  their  respective  velocity  profiles. 

3.3  Effect  of  Free  Stream  Reactant  Concentration  Ratio 
3.3.1  Product  Formation 

The  concentration  ratio,  €,  is  defined  in  equation 
A4  as  the  ratio  of  the  free  stream  osone  concentration 
to  nitric  oxide  concentration.  In  this  series  of  runs 
the  nitric  oxide  concentration  was  held  fixed  at  0.5% 
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while  the  ozone  concentretlon  verled  from  0.5  to 
7Z,  giving  e  range  of  C  from  1  to  14.  Free  stream 
velocities  were  5  end  25  ms'1  and  the  moveable  side 
walls  were  adjusted  for  zero  pressure  gradient  with 
zero  reaction.  The  small  negative  pressure  gradient 
(i.e.,  free  stream  acceleration)  caused  by  the  volisne 
expansion  in  the  reacting  case  was  negligible  at  these 
temperature  rises.  The  measured  variable  was  maximum 
mean  temperature  rise  obtained  in  a  traverse  100  am 
downstream  of  the  splitter  plate.  This  was  measured 
with  a  tungsten  cold  wire  (section  2.5.2),  which  was 
calibrated  during  the  run  with  a  thermocouple  at  a 
known  fixed  location  near  the  dividing  streaoiline.  A 
reference  junction  was  placed  in  one  free  stream.  Re¬ 
sponse  time  of  these  thermocouples  was  approximately 
0.4  sec  in  air  at  5  ms'1.  (The  faster  thermocouples 
described  in  2.5.2  were  not  available  at  this  stage 
of  the  experiment). 

The  results  obtained  are  plotted  in  Figure  3.15 
and  are  compared  with  the  predicted  values  obtained 
from  the  non-reacting  PDFs  measured  by  Konrad,  as  des¬ 
cribed  in  Appendix  A.  It  is  apparent  that  the  maximum 
mean  temperature  rise  predicted  from  the  PDF  is  very  close 
to  the  measured  value  at  these  low  temperature  rises. 

It  is  also  evident  that  the  temperature  rise  does  not 
appear  to  asymptotically  approach  a  limiting  value,  at 
least  within  this  range  of  £• 


The  existence  of  an  asymptote  can  be  demonstrated 
matheoiatlcally  in  the  case  of  infinite  reaction  speed. 
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Given  the  mixture  fraction  PDF,  Pc(c),  and  the 
temperature  rise  as  a  function  of  mixture  fraction  T(c) 
the  mean  temperature  rise  at  a  point  is  given  by 

AT  -  {  Pc(c)T(c)dc  (3.6) 

J0 

where 

AHf 

T(c)  "  c  C0  NO  -c7  f°r  C  -  C8 

(3.7) 

°r  AH. 

(1  -  c)  C  -7—  for  c  >c 
0  03  Cp  8 

c  is  the  stoichiometric  mixture  fraction,  and  is  shown 

8 

to  be 

C 

cs  ■  c - Vc -  "  T+J  **  Appendix  A. 

0  NO  0  03 

Hence 


— 

C 

0  NO  cp 


(c)dc  +  C 


(1  -  c)  p  (c)dc 
c 


(3.8) 


Ue  require  the  limit  as  5  m  of  AT*,  i.e.,  the  limit  as 
■  1.  Replacing  with  1  -  e  we  wTite 


_  1-c  .  _  1 

11a  t*  •  11m  |f  c  p  (c)dc  +  — ~  f  1  -  c)p  (c)dc 

LJ0  c  c  J 

0  l-e 


(3.9) 
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Given  that  p(c)  is  a  probability  density  function 
whose  integral  is  therefore  unity,  the  largest 
value  of  AT*  would  occur  in  the  limit  as  5  +  ■, 
only  if  pfi(c)  were  a  delta  function  arbitrarily 
close  to  c  ■  1.  i.e. ,  supposing  Pc(c)  ■  <5(1  -  a) 
then  AT*  ■  ^  for  c  >  a  and  1  -  a  for  e  >  a. 

Hence,  the  maximum  possible  value  of  AT*  is  one 
which  corresponds  physically  to  a  perfectly  mixed 
stoichiometric  mixture  in  which  €  approaches  infinity. 

In  Appendix  D,  it  is  shown  that  the  chemical  re¬ 
action  time  is  always  decreasing  with  Increasing  re¬ 
actant  concentrations.  This  leads  to  the  possibility 
of  a  Damkohler  number  effect,  since  at  the  lower  ozone 
concentrations  used  in  this  experiment  the  reaction 
time  becomes  comparable  with  the  larger  turbulent  time 
scales.  The  PDF  solution  assumes  of  course  an  Infinite 
reaction  speed  and  the  effect  of  finite  Damkohler 
number  would  be  to  give  a  lower  temperature  than  pre¬ 
dicted  by  the  PDF  solution.  This  may  account  for  the 
average  temperature  being  somewhat  below  the  predicted 
values  from  the  PDF  at  low  values  of  £  in  Figure  3.15. 
3.3.2  Skew  of  the  Mean  Temperature  Rise  Profile 

Figure  3.16  plots  three  temperature  profiles  obtained 
over  a  wide  range  of  concentratlonratio,  showing  a  skew 
in  the  profile  toward  the  minority  species  side.  For 
the  sake  of  comparison  individual  profiles  have  been 
scaled  to  the  same  height.  Figure  3.17  shows  the  same 
profiles  inferred  from  the  non-reacting  PDF  solutions 
with  the  same  scaling.  A  detailed  comparison  between 
experimental  and  predicted  profiles  is  made  in  Section  3.2 


Temperature  Rise  (Arbitrary  Scale) 
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but  this  discussion  concerns  the  effects  of  changing 
£.  The  three  arrows  at  the  top  of  the  figure  indicate 
the  position  of  the  peak  of  the  experimental  profiles 
in  3.16.  The  rightward  shift  of  the  Inferred  profiles 
(discussed  in  Section  3.2)  is  evident  here,  but  the 
relative  displacement  due  to  changing  £  are  apparently 
very  similar.  The  magnitude  of  this  skew  as  a  function 
of  £  becomes  meaningful  when  it  is  compared  to  the 
skew  obtained  in  a  laminar  diffusion  shear  layer  (or  a 
corresponding  turbulent  shear  layer  if  the  molecular 
diffusivltles  are  replaced  by  assumed  analogous  eddy 
dif fusivities) . 

The  solution  of  the  non-reacting,  uniform  density 
laminar  diffusion  shear  layer  in  the  temporal  frame¬ 
work  becomes  the  solution  of  the  one  dimensional  species 
conservation  equation  with  diffusion 


which  has  the  similarity  solution 


(3.10) 


ni  erf  (n)  +  1  (3.11) 

n  2 


where  n  is  the  similarity  variable 


n  n 

Since  the  number  density  is  constant,  -*■  ■  l - . 

n  n 

Knowing  the  mixture  fraction  for  any  y  it  is  simple  to 
infer  the  temperature  rise,  using  equation  3.7,  and  arrive 


•t  •  temperature  profile  for  any  given  £»  as  shown 
in  Figure  3.18.  This  analysis  Ignores  the  effect 
of  heat  on  the  density  and  diffusivltles,  but  Is 
adequate  for  the  present  purpose.  In  order  to  compare 
these  temperature  profiles  with  the  experimental  ones, 
the  n  coordinate  was  scaled  by  the  maximum  slope  thick¬ 
ness  of  the  species  concentration  profile  (equation  3.11), 
In  the  same  way  as  the  experimental  profiles  were  scaled 
by  the  maximum  slope  thickness  (6^)  of  the  velocity 
profile. 

Figure  3. 19  shows  the  normalized  temperature  pro¬ 
files  given  by  the  diffusion  model  for  the  same  values 
of  £  as  the  three  cases  in  Figure  3.18  and  3.17.  The 
arrows  indicate  the  positions  of  the  peaks  of  the  ex¬ 
perimental  profiles  and  show  that  the  diffusion  model 
gives  a  profile  skew  about  twice  as  large  as  the  skew 
actually  measured  in  the  turbulent  flow. 

This  result  is  another  clear  demonstration  of  the 
role  played  by  the  actual  PDF  and  the  misleading  picture 
that  might  be  obtained  if  the  turbulent  flow  were 
thought  of  as  being  in  some  ways  analogous  to  a  laminar 
problem  but  with  eddy  diffusivltles  replacing  molecular 
diffusivltles. 

3.4  Effect  of  Heat  Release 

3.4.1  Velocity  Profile  and  Entrainment 

It  was  observed  from  the  photographs  that  the 


spreading  rate  of  the  layer  with  heat  release  did  not 
appear  to  be  greatly  changed.  The  velocity  profiles 
presented  in  Figure  3.20  give  an  example  of  the  apparent 
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Fig.  3.18  Normalized  Temp  Rise  Profiles. 
Laminar  Diffusion  Model. 
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slight  thinning  of  the  layer  at  a  particular  downstr 


location  with  Increasing  heat  release.  Figure  3.21 

plots  6u/(x  -  xQ)  against  maximum  mean  temperature 

rise  for  layers  of  various  concentration  ratios.  It 

is  not  strictly  correct  to  interpret  this  quantity  as 

the  overall  spreading  rate  because  the  value  of  xq, 

the  displacement  of  the  effective  origin  of  the  fully 

developed  layer  from  the  splitter  plate  edge,  is 

assumed  constant,  and  was  measured  only  in  the  non- 

reacting  flow.  Its  determination  in  each  reacting 

flow  would  have  required  a  number  of  identical  runs  for 

each  value  of  heat  release,  which  was  not  a  practical 

proposition  in  the  present  rig.  Clearly,  the  effective 

origin  will  depend  on  the  growth  history  of  the  laminar 

part  of  the  layer,  a  process  which  is  influenced  by 

heat  release  as  indicated  by  the  photographic  evidence. 

A  study  of  a  number  of  photographs  of  the  flow  with 

various  heat  release  indicated  that  the  change  in  xQ 

would  have  been  of  the  order  of  7  mm  at  8  *  0.5. 

max 

(8  _  ■  AT  /T  . )  This  corresponds  to  a  reduction 

max  max  • 

of  (x  -  xQ)  by  about  6%  over  the  range  of  these 
experiments,  which  implies  that  at  any  fixed  value  of 
x,  about  two  thirds  of  the  perceived  thinning  of  the 
layer  is  due  to  transition  delay,  and  one  third  due  to 
reduced  growth  rate. 

Even  though  the  growth  rate  is  only  slightly 
reduced  it  is  obvious  that  the  entrainment  of  mass  into 
the  layer  has  been  strongly  reduced,  as  shown  in  Figure  3.22 
which  plots  the  normalized  mass  flux  against  8  _  . 
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The  Integration  limits  for  the  mass  flux  were  chosen 
6 

y 

at  y  ■  i  y  ,  these  points  being  relatively  Insensitive 
to  local  purturbations  in  the  mean  velocity  or  temp¬ 
erature  profile.  Inspection  of  Figures  3.21  and  3.22 
indicates  that  at  0  .  *0.5,  corresponding  to  a  max- 
lrnum  mean  density  decrease  of  33%,  the  mass  flow  is 
reduced  by  approximately  10%  due  to  thinning  of  the 
layer,  and  a  further  15%  due  to  the  density  decrease, 
a  total  decrease  of  25%,  As  explained  above,  about 
two  thirds  of  the  thinning  appears  to  be  caused  by  the 
heat  release  delaying  the  early  developmental  stages 
of  the  turbulent  shear  layer. 

3.4.2  Mean  Temperature  Rise 

Although  it  is  clear  that  the  total  entrainment 
is  reduced  by  heat  release,  it  is  not  obvious  if  the 
PDFs  have  changed  or  not.  The  temperature  rise  at  a 
point  is  determined  entirely  by  the  mixture  fraction 
PDF  at  that  point.  It  is  easy  to  show,  from  the  analysis 
presented  in  Appendix  A,  that  if  the  PDF  is  unchanged, 
and  £  is  fixed,  the  mean  temperature  must  be  a  linear 
function  of  the  reactant  concentration.  Due  to  exp¬ 
erimental  constraints,  it  was  not  possible  to  make  a 
large  number  of  runs  with  precisely  the  same  £.  How¬ 
ever  it  is  not  necessary  to  be  confined  to  a  fixed 
location  and  fixed  £  to  observe  evidence  of  the  change 
in  the  PDF  with  Increasing  heat  release.  This  can 
be  seen  by  plotting  6  from  the  experimental  profiles 

ul&X 

against  that  from  the  PDF  inferred  profiles  (Figure  3.23) 
for  the  same  reactant  concentrations.  The  points  shown 


Fig.  3.2)  Comparison  of  Measured  0  with  0  inferred  from 

max  max 

Konrad's  Non-reacting  PDFs.  (Nitrogen  Shear  Layer) 
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correspond  Co  a  vide  range  of  C,  and  the  excellent 
agreement  up  to  6  ■  0.35  tends  to  confirm  the 

Hill 

validity  of  the  PDFs  near  the  position  of  6 

ux 

Beyond  the  point  where  0 _  ■  0.35,  the  heat  re- 

max 

lease  Is  evidently  no  longer  passive  -  the  ectual 
PDF  begins  to  alter  significantly  and  the  non¬ 
reacting  PDF  is  no  longer  successful  in  predicting 
the  temperature  rise.  In  principle,  the  procedure 
used  in  Section  3.2  to  determine  a  simplified  PDF 
from  a  series  of  temperature  measurements  is  also 
applicable  to  solve  for  the  PDF  in  the  flow  with 
high  heat  release.  However,  this  would  require 
taking  a  series  of  runs  all  with  the  same  high  heat 
release  and  a  wide  range  of  5,  a  task  outside  the 
capabilities  of  the  present  apparatus,  and  indeed 
probably  not  feasible  with  the  ozone-nitric  oxide 
reaction  due  to  the  instability  of  concentrated  ozone. 
Broadly  speaking,  two  parameters  which  characterise 
any  PDF  are  mixedness  and  mean  mixture  fraction.  If 
we  assune  that  Che  entrainment  ratio  is  not  significantly 
changed  by  the  heat  release  (hence  maintaining  the  same 
mean  mixture  fraction)  then  the  fact  that  the  temperature 
does  not  achieve  the  value  predicted  by  the  non-reacting 
PDF  implies  that  the  molecular  scale  mixedness  must  have 

decreased  when  8  >  0.35.  We  note  that  this  is  con- 

max 

slstent  with  the  flow  visualisation  presented  in  Section 
3.1,  where  it  was  observed  that  the  fine  three  dimensional 
scales  were  progressively  suppressed  by  the  heat  release. 
Should  this  be  the  major  cause  of  the  deviation  shown 
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In  Figure  3.23,  the  crucial  rola  of  the  small  scales  in 
mixing  is  apparent.  It  is  interesting  to  compare  the 
result  of  Konrad,  showing  that  the  mlxednese  at  Reynolds 
numbers  above  the  critical  value  (l.e.  after  transition) 
is  about  25Z  higher  than  it  was  prior  to  transition.  The 
increase  in  mixedness  coincided  with  the  appearance  of  the 
small  scales. 

3.5  Effect  of  Free  Stream  Density  Difference  on  Product 

Formation 

It  has  been  established  (Brown,  1974)  that  the  relative 
entrainment  rates  of  the  free  stream  fluids  in  a  turbulent 
mixing  layer  is  dependent  upon  their  relative  densities.  Ir¬ 
respective  of  whether  it  is  on  the  high  velocity  or  low  veloc¬ 
ity  side,  more  of  the  lower  density  stream  is  entrained  than 
of  the  higher  density  stream.  This  property  is  clearly  evident 
in  the  PDF  of  Konrad  (Figure  A2)  in  which  the  most  probable  con¬ 
centration  at  almost  all  points  across  the  layer  is  skewed  to 
the  helium  side  (c.f.  the  constant  density  PDF  of  Figure  Al). 

One  way  of  looking  at  the  effect  is  to  note  that  the  entrain¬ 
ment  rate  of  helium  is  about  3  times  that  of  nitrogen,  in  con¬ 
trast  to  the  uniform  density  layer  where  the  ratio  of  entrain¬ 
ment  rate  is  about  1.  The  question  as  to  why  the  variable 
density  has  such  a  significant  effect  on  the  PDF  is  considered 
later  in  this  section.  Vhether  this  effect  will  lead  to  an 
increased  or  decreased  product  concentration  (analogous  to 
temperature  rise)  depends  on  the  concentration  ratio,  €,  of  NO 

and  0  in  their  respective  free  streams,  that  is,  on  the 

3 

concentration  c#  at  which  NO  and  0^  are  in  chemical  equivalence. 
Generally  speaking,  the  highest  mean  product  concentration  is 
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expected  when  the  ratio  of  NO  and  0  free  stream 

3 

concentrations  corresponds  to  the  entrainment  ratio 
of  the  free  streams.  This  Is  because  at  the  most 
probable  concentration  (which  corresponds  roughly 
with  the  entrainment  ratio) ,  the  chemical  equiva¬ 
lence  will  be  near  unity. 

To  demonstrate  the  crucial  effect  of  a  free 

stream  density  ratio  change,  two  runs  were  made  with 

the  same  reactant  concentrations,  free  stream  velocities 

and  streamwlse  pressure  gradient  (zero).  The  only 

difference  was  to  change  the  high  speed  stream  from 

nitrogen  to  helium.  Nitrogen  was  used  for  the  low 

speed  stream  in  both  cases,  and  the  reactant  concentration 

ratio  £  was  1.  Despite  an  increase  In  the  flame  sheet 

temperature  by  a  factor  of  1.2+  in  the  helium-nitrogen 

case,  due  to  the  reduced  value  of  C  ,  the  maximum  mean 

P 

temperature  rise  in  the  turbulent  layer  fell  to  0.6  of 
that  In  the  nitrogen-nitrogen  case.  Figure  3.24  gives 
the  experimental  profiles.  This  result  Is  in  good 
agreement  with  the  Inferred  temperatures  calculated 
from  Konrad's  PDFs  (see  Section  3.2.1).  The  mean 
mixture  at  the  dividing  streamline  now 


t  If  the  Lewis  number  D/a  is  one  and  the  diffusivlties 
of  all  species  are  equal. 
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strongly  favours  helium  (c  -  0.8,  according  to  Figure 
3,  Konrad  1976)  because  of  the  Increased  entrainment  of 
the  lower  density  gas.  The  position  of  the  mixture 
leading  to  the  maximum  temperature  (about  1:1  since  £ 

Is  1)  Is  now  moved  strongly  toward  the  nitrogen  free 
stream,  causing  a  shift  to  the  left  of  the  He/N 

2 

profile  on  Figure  3.24. 

The  result  is  in  direct  contrast  to  that  which  would 
be  obtained  In  a  laminar  diffusion  (or  the  analogous 
eddy  dlffuslvlty)  model  of  the  shear  layer,  which  would 
predict  a  higher  temperature  in  the  He/N^  case.  The 
rigorous  solution  for  the  laminar  diffusion  model  with 
unequal  densities  is  discussed  In  Brown  and  Wallace  (1981) 

In  Section  3.1.1  the  significance  of  Interacting 
density  gradient  and  acceleration  fields  In  causing  a 
redistribution  of  the  primary  vortlclty  by  barocllnlc 
torque  action  was  noted.  A  similar  analysis  of  this  non- 
uniform  density  layer  provides  an  explanation  for  the 
asymmetric  appearance  and  entrainment  rate.  The  dotted 
line  in  Figure  3.25  represents  an  isopycnic  surface 
between  the  helium  and  nitrogen  free  streams,  which  Is 
normal  to  the  density  gradient,  (Patnalk,  Sherman  and 
Corcos,  1976).  The  barocllnlc  torque  term  (refer  to 
Section  3.1.1)  on  this  line  is  seen  to  be  strengthening 
the  general  vortlclty  in  the  region  marked  A  and  weaken¬ 
ing  it  at  B.  The  Induced  velocity  (arrow  C)  in  the 
helium  free  stream  is  consequently  greater  than  that 
In  the  nitrogen,  with  the  result  that  a  larger  volvne 


of  helium  enters  the  structure,  thereby  distorting  it 
to  the  shape  observed  in  Figure  3.2. 

3.6.'  Mean  Temperature  Rise  in  Argon  Shear  Lavers 

Obviously,  carrying  the  reactants  in  a  gas 
with  a  lower  specific  heat  should  produce  a  higher 
temperature  rise  if  all  other  flow  parameters  are  the 
same.  For  this  reason  it  was  proposed  to  use  argon 
instead  of  nitrogen  as  the  carrier  gas  in  order  to 
achieve  a  higher  temperature  rise  for  a  given  reactant 
concentration,  thereby  enabling  the  heat  release  effect 
shown  in  Figure  3.23  to  be  checked  to  higher  temperature. 
Upon  doing  the  experiment,  however,  it  was  most  surprising 
to  find  that  this  expected  increase  in  temperature  (an 
extra  402)  did  not  occur.  Indeed,  performing  two  con¬ 
secutive,  identical  runs,  except  for  merely  swapping 
nitrogen  for  argon  (a  simple  procedure  on  this  facility), 
produced  a  slightly  lower  maximum  mean  temperature  in 
the  argon  layer  than  in  nitrogen!  Numerous  checks  of 
the  facility  and  probes  were  made,  and  a  large  number  of 
runs  taken.  The  results  are  shown  in  Figure  3.26,  which 
is  directly  comparable  with  Figure  3.23.  The  velocities 

were  5  and  25  ms*1  and  pressure  1  atm.  0  was  measured 

max 

at  x  ■  100  mm. 

At  the  time  of  writing,  it  is  still  unclear  why  the 
temperature  rise  in  the  argon  shear  layer  falls  so  short 
of  that  predicted  by  the  non-reacting  PDFs  of  Konrad. 

The  most  likely  explanation^  appears  to  be  one  connected 

t  First  suggested  to  me  by  Dr.  Keith  King,  University  of 
Adelaide. 


with  the  chemical  kinetics.  It  is  reasonably  wall  established 

that  the  oxone-nltrlc  oxide  reaction  is  basically  a  simple, 

one  step,  blaolecular  reaction  (Johnson  &  Crosby,  1954),  but 

the  appearance  of  the  energy  released  (as  a  temperature)  can 

* 

be  speed  limited  by  the  rate  at  which  the  activated  NOa  trans¬ 
fers  its  energy  no  the  carrier  gas.  l.e.  the  relevant  reactions 
so  far  as  the  temperature  rise  is  concerned  are 

1.  NO  0,*Nc£  ♦  0* 

2.  NO*  +  M  *  NO,  +  M  +  heat 

The  reaction  time  quoted  in  calculating  the  Damkohler  number 
(appendix  D)  refers  only  to  the  first  reaction,  because  it  was 
determined  by  measuring  the  disappearance  of  ozone.  The 
rate  of  reaction  (2)  has  been  assumed  to  be  fast  compared  to 
reaction  (1)  making  (1)  the  rate  limiting  step  on  the  molecular 

scale.  The  rate  of  reaction  (2)  depends  on  the  collision 

* 

kinetics  of  the  molecules  N0a  and  M,  and  in  particular  on 
whether  M  is  bimolecular  (e.g.  nitrogen)  or  atomic  (e.g. 
argon),  i.e.  on  the  cuopling  efficiency.  This  coupling  cannot 
be  predicted  by  the  cross  sectional  analysis  model  of  molecular 
interaction.  It  is  conceivable,  (although  no  direct  confirmation 
has  been  found  in  the  literature  to  date)  chat  the  coupling 

1r 

of  the  vibratlonally  excited  modes  of  N02  to  the  molecule 
M  will  be  more  efficient  if  N  is  nitrogen  than  argon. 

If  this  is  the  case,  then  the  observed  discrepancy 
is  the  result  of  a  Damkohler  number  effect. 
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4.  CONCLUDING  REMARKS 

The  application  of  the  unique  flow  facility 
described  has  begun  to  shed  sone  light  on  the  complex 
Interaction  between  chemical  heat  release  and  turbulent 
mixing,  perhaps  the  least  well  understood  aspect  of 
combustion  at  this  point  In  time. 

It  has  been  demonstated  that  although  the  heat 
release  has  produced  a  volume  expansion  In  the  layer, 
the  entrainment  rate  of  free  stream  fluid  has  been 
reduced  by  about  the  same  amount,  causing  the  overall 
growth  rate  to  remain  substantially  unchanged.  Further 
understanding  may  determine  if  this  is  coincidental 
or  has  some  mechanistic  explanation.  According  to  the 
pictures,  the  mixedness  of  the  fluid  in  the  layer  is 
apparently  reduced  by  the  heat  release  as  the  small 
scales  are  progressively  reduced.  This  is  confirmed 
by  mean  temperature  measurements  which  show  that  the 
relationship  between  mean  temperature  rise  and  reactant 
concentration  (for  a  fixed  reactant  concentration  ratio) 
is  not  linear  beyond  a  temperature  rise  of  IOO°K.  Experi¬ 
ments  ranging  to  much  higher  temperature  rises  are  needed 
to  fully  observe  this  trend  and  elucidate  its  implica¬ 
tions  on  the  shape  of  the  mixture  fraction  PDFs. 


APPENDIX  A 


Numerical  Prediction  of  Temperature  Rise  Using  Mixture 
Fraction  PDFs 

A1  Konrad's  PDF  Data 

The  extent  of  molecular  scale  mixing  in  the  shear 
layer  was  measured  with  a  Brown-Rebollo  (1972)  con¬ 
centration  probe  having  a  spaclal  resolution  of  200  pm 
and  a  frequency  response  to  65  kHz.  Konrad  obtained 
mixture  fraction  PDFs  at  seven  positions  across  the 
layer  for  both  the  uniform  and  non-uniform  density 
case.  These  are  shown  in  Figures  A1  and  A2. 

It  is  possible  to  use  this  data  to  predict  the 
average  and  fluctuating  species  concentrations  in  the 
flow,  and  in  particular  find  the  product  produced 
in  a  reacting  flow.  To  do  this  we  need  to  know,  for 
a  perfectly  mixed  particle,  the  resultant  product 
concentration  as  a  function  of  the  mixture  fraction 
and  the  two  free  stream  reactant  concentrations.  This 
is  derived  in  Section  A2  below.  The  procedure  to 
calculate  the  mean  product  concentration  and  temperature 
rise  is  given  in  A3. 

Such  a  prediction  involves  the  following 
assumptions: 

(1)  No  effect  of  heat  release  on  the  fluid  dynamics. 

(2)  An  infinite  Damkohler  number,  i.e. ,  mixing 
limited  reaction. 

Assumption  (1)  is  reasonably  supported  up  to  reactant 


00 

AlligVQOdd 


Fig.  A2  Probability  Density  Function;  r  »  0.38 
(From  Konrad). 
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concentrations  of  about  0.51  (only  51  change  in  density) .  Clearly 
this  is  not  satisfied  in  the  reacting  flows  but  by  comparing  pre¬ 
dictions  based  on  this  assumption  with  the  actual  results  it  is 
possible  to  determine  the  extent  to  which  isothermal  mixing  is 
affected  by  the  exothermic  reaction.  Refer  Figure  3.23.  The 
Damkohler  number  is  estimated  in  Appendix  D,  where  it  is  shown  that 
the  reaction  rate  is  limited  by  the  mixing  process  in  this 
experiment  except  where  the  reactant  concentrations  are  less  than 
1%.  Most  of  the  runs  in  this  work  were  mixing  limited. 

A2  Perfect  Mixing  Analysis 

The  dominant  reaction  occuring  in  the  flow  is 
0,  +  NO  *  NO*  +0* 

Consider  two  volumes  of  gas  each  with  a  total  number  of  nt  and  n* 

molecules  respectively.  The  nt  molecules  Include  nu_  molecules 

NO 

of  NO  in  otherwise  inert  molecules,  and  likewise  the  n*  molecules 
include  n^  molecules  of  ozone.  We  define  the  free  stream  con¬ 
centrations  of  NO  and  Oj  as 


C, 


03 


(A2) 


and  a  mixture  fraction 

c 


U| 


ni  ♦  n* 


(A3) 


We  also  define  the  free  stream  reactant  concentration  ratio 


(A4) 


5 


Now  the  two  volumes  of  gas  are  mixed  completely, 
i.e.,  the  reaction  proceeds  until  one  of  the  re¬ 
actants  is  completely  exhausted.  The  concentration 
of  the  product  NO  is  given  by 


or 


1  2 


_ _ n™- 

w02  n  +  n 
1  2 


which  become 


if  >  n 


"NO  *  "03 


if  n«,  > 


03  >  "NO 


(A5) 


'N02  -  c 


So 


or 


if  nQ3  >  nNQ 


-N02  -  (1  -  c)C03  if  nN()  > 


*03 


(A6) 


At  the  point  where  nQ3  •  n^g,  the  reactants  are  present 
in  their  stoichiometric  ratio,  and  the  stoichiometric 
mixture  fraction  is  given  by 

•  (1  -  c.)C03 


i.e. 


and 


'03 


8  C03  +  CN0 


's  C  +  1 


(A7) 


(A8) 


This  mixture  fraction  leads  to  the  maximum  achievable 
product  concentration 


CN02nax  •  CsSo 


(A9) 


So2  as  a  function  of  c  is  plotted  in  Figure  A3.  In 


perfect  mixing  then,  the  operetlng  point  aovee 
along  the  path  JOM,  achieving  a  maximum  CNQ2  at 
a  mixture  ratio  c#  defined  by  equation  A8. 

Hence 

S02(c)  "  So(c  H(cs  *  c>  ♦  *<l  -  C)H(C  ~  *.»  (A10) 

where  H(x)  i a  the  Heaviside  step  function. 

In  the  case  where  both  free  streams  have  the 
sane  specific  heat,  temperature  rise  becomes  the 
direct  analogue  of  product  concentration  since 


T(c)  • 


"K02AHf 

(n  +  n  )C 
1  2  P 


C~  “  COMt*  Sl02(c) 


(All) 


where  AHf  is  the  heat  release  per  mole  of  reaction 

and  C  is  the  average  specific  heat  per  mole  of  the 
P 

mixture  at  constant  pressure.  The  abscissa  of  Figure 
A3  could  equally  well  be  temperature  rise. 

In  the  case  where  the  two  incident  streams  have 
different  specific  heats,  the  specific  heat  of  the 
mixture  is  a  function  of  mixture  ratio.  Hence 


T(c) 


tHICH02(c> 

C  (c) 

P 


(A12) 


where 
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A3  Calculation  of  Product  Concentration;  and  Temperature 
Riga  Profllaa 

The  seen  product  concentration  la  calculated  froa 
lta  own  PDF  ualng  the  Identity 


■f 


^02  PN02(CN02)dCN02 


°N02 


(A14) 


•  function  of  c  (equation  (A10)  which  Is 
now  made  the  variable  of  Integration 


_  f  .  r  {r  \  dCW02 

Sl02  "  J  "02(c)  PN02'S02;  dc 
c  , 


dc 


(A15) 


The  quantity  in  square  brackets  is  the  mixture  fraction 


PDF  Pc(c),  by  Identity  (e.g.  Fry  1965). 
Hence 


’«('>dc  <A16> 

Similarly, 

*  f  T(c>  pc<c>  dc  (A17) 

Jo 

where  T(c)  is  defined  In  equations  All  and  A12.  To 

construct  a  temperature  profile,  equation  A17  was 

solved  with  each  of  the  seven  PDFs  P  (c)  measured  by 

c 

Konrad.  This  was  done  numerically  with  a  simple 


FORTRAN  program. 


APPENDIX  B 


B1  Calculation  of  the  Dividing  Streamline 

The  position  of  the  dividing  streamline  vas 
found  from  the  mean  temperature  and  velocity  profiles 
using  the  following  expression 

m 

|  (p  u  0^  -  0  «?)dy  ■  |  (p?*pu  U^dy 

y*  — 

This  estimation  is  necessarily  an  approximation  since 
a  number  of  higher  order  terms,  as  well  as  correlations, 
have  been  omitted.  The  expression  is  simply  obtained 
from  the  conservation  of  continuity  and  momentum  in  a 
strip  extending  to  infinity  each  side  of  the  dividing 
streamline. 
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APPBTOIX  C 


Digital  Filter  Characteristic* 

Tha  digital  filtar  was  uaad  to  smooth  both 
tha  temperature  and  valocity  signal  and  la  Identical 
to  that  used  by  Thomas  (1977)  previously  of  this 
department ,  who  provided  much  of  the  Information 
appearing  in  this  section.  It  was  described  by  Linn 
(1977). 

Physically,  the  filter  steps  through  the  data 
record,  averaging  a  fixed  span  of  samples  and  placing 
the  result  at  the  span  centre.  Z  plane  analysis  is 
used  to  obtain  the  frequency  response  and  phase  shift 
characteristic.  The  Z  transform  in  this  case  is 

1  -  xk  n 

Y(z)/X(z )  -  (j  :  |-) 

idie  re 

k  »  order  of  the  filter 
n  *  number  of  applications 

Y,X  ■  output  and  input  respectively 
,  JwP 

P  -  sample  period 
u  ■  circular  frequency 

which  represents  k  zeroes  equispaced  around  the  unit 
circle,  with  one  (at  u  ■  0)  being  cancelled  by  a 
coincident  pole.  The  time  domain  operation,  obtained 
by  inspection  from  equation  Cl  is  for  r  »  1. 


(Cl) 
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Yn  “  +  l  "  X-  *  X 


n  n  +  k 


l.«. 


Y  ■  Y  -  X  +X 
n  +  1  n  n  n+k 


(C2) 


The  record  hes  experienced  e  shift  of  (k  -  l)/2 

sample  periods  by  this  process,  or  after  n  applications 

of  the  filter,  a  total  shift  of  n(k  -  l)/2  sample 

periods.  By  choosing  an  even  number  for  n,  a  shift 

of  an  Integral  nianber  of  sample  periods  is  assured. 

The  shift  is  obviously  Independent  of  frequency,  and 

is  readily  nulled  by  an  equal  shift  of  opposite  sign 

after  the  filtering  process.  Choosing  a  higher  value 

of  n  also  has  the  advantage  of  giving  the  filter  a 

shorter  cutoff.  The  filter  is  recursive  because  it 

is  using  a  previous  output  (Yr)  to  derive  the  new 

output  (Y  +1).  This  feature  causes  the  RH  side  of 
n 

equation  C2  to  have  only  3  terms  and  permits  rapid 
computation. 

Although  the  actual  filter  cutoff  is  only  of 
secondary  consideration  for  present  purposes ,  a  typical 
frequency  response  is  plotted  in  Figure  Cl.  The 
computations  of  the  curve  was  extended  to  a  frequency 
of  5/2P  to  graphically  show  the  Nyqulst  folding  phenomenon 
within  the  filter  at  1/P,  2/P.*  Such  folding  is  a  danger 
since  it  causes  high  frequency  energy  to  appear  as  low 
frequency  energy  after  the  digital  processing.  The  only 
way  of  avoiding  it  is  .to  ensure  (by  using  analogue  fitters 


t  A  good  discussion  of  this  phenomenon  is  to  be  found  in 
Ben da t  &  Piersol  (1971). 


[I 


Frequency  Hz 


Fig.  Cl  Relative  Frequency  Response  of  Low  Pass  Filters 

(a)  n  «  2,  k  -  2 

(b)  n  ■  2,  k  ■  10 
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if  necessary)  that  the  analogue  signal  contains, 
neglible  energy  content  above  frequency  of  1/2P. 

The  sampling  periods  used  were  selected  on  this 
basis  and  checked  by  Inspection  of  an  expanded  plat 
of  the  dlgital-to-analogue  converted  signal  to  en¬ 
sure  that  an  acceptably  "clean"  analogue  signal 
could  be  obtained  by  connecting  all  the  plotted 
points  with  a  straight  line. 
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APPENDIX  D 


Calculation  of  the  Damkohler  Number 

The  purpose  of  this  section  Is  to  calculate  a 
relative  time  scale  for  the  chemical  reaction  and  com¬ 
pare  It  to  the  range  of  mixing  time  scales  In  the  flow. 
Mixing  In  the  shear  layer  occurs  over  a  wide  range  of 
length  and  time  scales,  the  largest  being  the  scale  of 
the  width  of  the  shear  layer  (because  no  larger  structure 
could  possibly  be  Incorporated)  and  the  smallest  being 
the  scale  at  which  molecular  diffuslvlty  becomes  dominant 
enough  to  smear  the  inhomogenities  into  molecular 
homogenlety. 

D1  Mixing  Time  Scales 

The  largest  time  scale,  called  the  convection  time 

scale  t  Is  derived  from  the  mean  velocity  (U  +  U  )/2  and 
c  12 

the  layer  width  at  the  cross  section  of  Interest. 

26 

T  -  - ±- 

C  w  +  U 
l  2 

In  this  experiment,  at  x  ■  100  mm, 

6  ■  20  mm 

u  • 

U  ■  25  m/s 
1 

U  »  5  m/s 

2 

hence  the  typical  convection  time  is 
Tc  ’  L'Vy2 "  660  *  l0~6  #ec 
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The  snalleat  turbulent  length  scale  is  the 
Kolmogorov  microscale 

X  -  (v2/e)k 


v  -  kinetic  viscosity 
e  ■  K.E.  dissipation  rate  per  unit  mass. 
The  dissipation  rate  can  be  taken  as 

e  -  (U  -  0  )3M 
1  2  “ 


If  the  molecular  dlffusivlty  is  D,  then  the  time 
scale  appropriate  to  the  smearing  of  these  smallest 
eddies  is 

X2 


which  becomes 


Re**  *  ~ 


u 

D  AU 


The  Schmidt  no.  ^  for  the  species  in  this  flow  is  about 
unity.  Hence 


For  most  of  the  flow  in  this  case 
A  U  6 


Re  ■ 


20  x  .02 
1.6  x  10~5 


2.5  x  10“ 


TD  ■  -  6.3  x  10  6  secs. 

D2  Reaction  Speed  Time  Scale 

The  reaction  rate  of  the  second  order  reaction  is 


d  C 
“dt 


01 


•  k  C 


NO 


C 

03 
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The  ratio  of  mixing  time  to  chemical  reaction  time  is 

Tmix 

known  as  the  Damkohler  number,  a  ■  - -  .  The  mixing 

chem 

rate  is  limited  by  the  slowest  of  the  steps,  i.e.  t£. 
The  value  of  a  is  therefore  about  16  at  5%  reactant 
concentration,  falling  to  1  at  0.3%  reactant  concentra¬ 
tion.  Thus  the  rate  of  product  formation  is  controlled 
by  the  mixing  process  (the  large  scale  time)  unless  the 
reactant  concentrations  are  somewhat  less  than  1%. 


Orone-Nitric  Oxide  Chemist 


The  following  ractions  are  known  to  occur  in 
the  range  300  to  1000  K .  Reaction  1  la  dominant 
between  300  and  500°K. 

Reaction  rate  kf  •  Mp*  exp  (E./RT) 


Reaction 

A 

cm3 

molecule-1 

see”1 

NO  +  0  -*■  NO  +0 

9.0  x  10  13 

3  2  2 

NO  +  NO  +0  *  NO  +  NO 

1.0  x  10-12 

2  2  2 

0  +  M  0  +  0  +  M 

6.4  x  10-10 

3  2 

0  +  NO  -►NO  +  0 

9.1  x  10-12 

2  2 

0  +  NO  -►NO  +  0  +0 

9.8  x  10-12 

3  2  2  2 

0  +  0  -►  0  +0 

1.5  x  10-11 

3  2  2 

N0  +  0  +  M*N0  +  M 

1.0  x  10-31 

2 

Note  also  for  reaction  1  AH^  *  -  47.5  kcal  mole-1 

•c  -  [NO  ] [0  ]/[N0] [0  ]'  5  x  1034 
2  2  3 
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